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METHOD FOR LASER INDUCED ISOTOPE ENRICHMENT 



Cross Reference to Related Applications 

[0001] This application is a continuation in part of U.S. application serial no. 
09/914,401 which is a 35 USC 371 filing of International Application No. PCT/USOO/06077 
filed March 9, 2000. 

Government Rights 

[0002] This invention was made with government support provided by the US 
Department of Energy through the Office of Basic Energy Science (Chemical Sciences) 
under contract no. DE-FG02-01ER and through the National Science Foundation under 
Contract Number PHY-8920108. The government has certain rights in the invention. 

Field of the Invention 

[0003] This invention relates to methods for isotope enrichment. 

Background of the Invention 
[0004] Selected isotopes are useful in a variety of applications. For example, it is 
often desired to have an isotopic marker present in a molecule to trace such molecule in 
biologic processes. Radioactive isotopes are used for curative purposes in cancer and for 
medical radiology. Isotopic enrichment is used by industry to control physical properties of 
materials such as thermal conductivity. Separation of isotopes is conventionally 
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accomplished by a variety of means relying on difference in mass of isotopes. One example 
of this utilizes a gaseous diffusion process to separate heavier components of an isotopic 
mixture from the lighter components by reason of their differing mass. Physical separation 
by fractionation to isolate a particular desired isotope is also known. Conventionally, the 
most popular procedure for isotope enrichment is the gaseous-diffusion method, which 
depends on the relative rates of gaseous diffusion through molecular barriers by different 
isotopes of a given element. Such processes typically operate at reduced pressure and 
elevated temperature. Often the difference between the mass of the isotopes is very small 
and although the lighter weight isotope diffuses more rapidly through porous barriers, the 
separation by each barrier is very, very small, therefore, many or a cascade of barriers are 
required. 

[0005] More recently, there has been an interest in attempting to cause isotope 
separation by multiple photon absorption and photoionization. These processes are labor 
intensive and difficult to accomplish. Typically, several iterations are required to achieve the 
high concentration of the desired isotopic component in the separated product. An example 
of separation by multiple photon absorption can be found in U.S. Patent No. 4,655,890 
where a mixture of molecules is irradiated with infrared laser radiation having a first 
wavelength which lies in close coincidence with the transition from the ground rotation- 
vibration level of the molecules which contain the isotopic variant of the element which is to 
be separated, and the next highest rotation-vibration level having vibrational excitation. This 
radiation selectively induces multiple photon absorption by those molecules. Then it is 
possible to cause the selectively excited molecules to undergo chemical reaction while other 
isotopic variance of the element remains substantially in their lower rotation-vibration levels 
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and consequently unreactive, thereby causing a product of the chemical reaction of the 
selectively excited molecules to be produced and separable from the unreacted molecules. 

[0006] In the parent of the present application, a novel process and system 
apparatus were described for separating isotopes of an element and causing enrichment of a 
desired isotope of an element in a material. In one aspect, a laser was used to modify or 
fabricate a material so as to produce a desired isotopic content in the material, which differs 
from that which naturally occurs. 

[0007] Isotope enrichment in laser ablation plumes is predominantly an ionic 
process. Such a process involves the interaction of energetic ions produced in the ablation 
plasma with self generated magnetic fields in the same plasma. Details about the theory of 
this self generated magnetic field process and its effects on the motion of isotopic ions in an 
ablation plasma continue to be studied and worked out. However, it is recognized that by 
maximizing the ionic component of the plasma, the number density of isotopes that undergo 
the enrichment process can be increased. 

[0008] Referring to Figures 19A and 19B, it can be seen that the ratio of ions to 
neutrals in a single ablation plasma plume is relatively small. Typically that ratio is on the 
order of 1% or less. To achieve a practical method for using this phenomenon as a way to 
harvest enriched isotopes, it would be desired to increase this ionic component to much 
greater levels. 

[0009] A single laser ablation plume is the result of a number of relatively 
complex processes occurring during the absorption of a laser pulse by a material and 
subsequent expansion of the gaseous/plasma plume. One part of this process involves the 
initial production of a very large percentage of ions if not complete ion content in the vapor 
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phase plasma. In other words, the system starts out as a fully ionized electron-ion plasma. 
Because of the high initial density and ensuing collisions in at the early vapor/plasma stage, 
however, a great deal of electronic neutralization occurs as the plasma expands away from 
the surface of the ablating material. This highly neutralized plasma then moves through the 
region of the self-generated magnetic field and the ionic component undergoes the isotope 
enrichment. 

[0010] A plasma having a density which exceeds a certain level, known as the 
"critical density", begins to reflect optical radiation rather than absorb it. The absorption 
properties of a plasma normally increase with the density of free electrons in it. However, 
when the density reaches the critical density the plasma becomes nearly totally reflective. 
This limits the amount of internal energy that can be added to a plasma by a laser pulse. In 
most cases where an ablation plasma is formed by a laser this critical density is achieved in 
the early part of the laser pulse impinging on the surface of the material. 

Summary of the Invention 

[001 1 ] The present invention provides a novel process and system apparatus for 
separating isotopes of an element and causing enrichment of a desired isotope of an element 
in a material. In one aspect, the invention utilizes lasers to modify or fabricate a material so 
as to produce a desired isotopic content in the material, which differs from that which 
naturally occurs. The reference to a condition different from that which naturally occurs 
refers to both enhancement and depletion. 

[0012] In a preferred embodiment, the present invention relies upon laser induced 
formation of a plasma by means of a pulsed laser beam. The material containing the 

4 



Attorney Docket No. 2115D-002245 

isotopes desired to be separated is placed in the beam path and is called a target. A plasma is 
formed from the target material by focusing an intense, short duration optical pulse from the 
laser onto the target. The concentrated energy contained in the focused and/or concentrated 
laser beam ionizes the target material, energizing its electrons and raising it to a temperature 
or otherwise exciting it to a condition whereby ions are produced, thereby generating a 
plasma. The plasma contains ions having varying isotopes, ion energy and charge state 
distributions. In one aspect, the ions are contained in a plasma plume generated when laser 
pulses are directed to a target. The process and system of the invention produce an unusual 
isotopic enhancement effect in the observed ion spectra. This effect is manifested as an 
enrichment of the Ughter isotope in zones of the ablation plume, as observed normal to the 
surface of the target material. 

[0013] In one aspect, the method comprises directing the laser beam to the target 
at an intensity and wavelength sufficient to generate the plasma comprising ionized isotopic 
species and to generate an internal electromagnetic scattering field within the plasma causing 
spatial separation of the ionized isotopic species. In another aspect, the ions are collected on 
a substrate in a manner which provides zones having isotopic distribution different from that 
which naturally occurs for a given element contained in the target. It is preferred that the 
pulses which generate the plume be produced by an ultrafast laser (femtosecond or 
picosecond). It is preferred that the laser generates the plume by ablation or laser induced 
breakdown (LIB) of the target material. 

[0014] In another aspect, the method comprises directing an initial pulse of the 
laser beam to the target at an intensity and wavelength sufficient to generate the plasma 
comprising ionized isotopic species and to generate an internal electromagnetic scattering 
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field within the plasma causing spatial separation of the ionized isotopic species. The dense 
plasma is then allowed to expand so that its volume increases and its electron density 
decreases. This brings the density of the plasma equal to or less than the critical density. 
During the expansion process a certain degree of neutralization occurs in the plasma which 
additionally reduces the free electron concentration. After allowing the expansion to occur, 
a time delayed second pulse of the laser beam is directed to the plasma to further spatially 
separate the ionized isotopic species. By correctly delaying the time of the second pulse, 
additional energy is absorbed by the plasma and the ionic component increases. 

[0015] The invention provides substantial advantages over conventional methods 
for separating isotopic components of an element. The invention further provides the ability 
to effectively enhance the isotopic content of a material as compared to the natural state. 
Products having desired isotope distribution are able to be produced from the process. The 
process and system are adaptable to commercial use and automated production. Therefore, 
isotope enriched products are obtainable by the methods and system of the invention. 

[0016] Objects, features, and advantages of the invention include, in addition to 
the foregoing, an improved method and system for separating isotopes of an element, and 
particularly for forming product material having a desired isotopic content different from 
that which naturally occurs. 

[0017] Another object is to provide a method for producing spatial isotopic 
separation. 

[0018] Another object is to provide a method and system which utilize lasers to 
achieve isotope separation or enrichment at certain ablation or scattering angles. 
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[0019] Another object is to provide a method and system which utilize lasers to 
achieve isotopic separation or isotopic enrichment of a product. 

[0020] Another object is to provide a method and system for producing products 
having a desired isotopic distribution. 

[0021] Another object is to provide a method and system which effectively and 
economically achieve isotopic separation and isotopic enhancement in a product material. 

[0022] Another object is to provide a method and system which applies the 
technique of isotopic separation to chemical element separation. 

[0023] These and other objects, features, and advantages will become apparent 
from the following description of the preferred embodiments and accompanying drawings. 

Brief Description of the Drawings 

[0024] Figure 1 is a schematic of a system for generating a plume containing 
ionized particles. It includes a deposition chamber which is preferentially under vacuum 
(but could be operated with a background gas at various pressures) and an analysis or 
collection chamber. An incoming laser beam is provided from an external laser source. 

[0025] Figure 2 is a graph showing the signal from the microchannel plate 
(MCP) of the ion energy analyzer in units of voltage as a function of arrival time after the 
laser pulse in microseconds. Figure 2 shows individual charge states and kinetic energies 
for the ionic isotopes of boron obtained from a boron nitride target in a vacuum. The data of 
Figure 2 was captured in the experimental apparatus as shown in Figure 1. The laser pulse 
was supplied by a Ti:Sapphire CPA system. 
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[0026] Figure 3 shows the charge state and energy distribution data for boron 
ions and its isotopes, as compiled from spectra such as that on Figure 2. Figure 3 has two 
parts. Figure 3A shows data obtained in the experimental apparatus of Figure 1 where the 
deposition chamber was under vacuum, and Figure 3B shows data obtained in the same 
experimental configuration and under the same conditions except that nitrogen gas was 
included in the chamber. In Figure 3 for the ion species of nitrogen, boron 10, and boron 
11, there is shown data obtained from a cross sectional area transverse to the direction of ion 
movement and normal to the surface of the target. The data shows ion density in ions per 
square centimeter as a function of ion kinetic energy expressed in keV. The data of Figures 
2 and 3 were obtained under the conditions of boron nitride ablation with laser pulses of 200 
femtoseconds, 780 nanometers, 4 millijoules (mJ), 100 micron spot size, corresponding to 
2.5 X 10^"^ watts per square centimeter. 

[0027] Figure 4 contains two plots labeled as 4 A and 4B. The data of Figure 4 
was obtained by ablation of boron nitride target at 6 nanoseconds pulse width, 1.06 micron 
wavelength, 60 millijoule energy per pulse, corresponding to an energy fluence of 760 J/cm^ 
based on 100 micron spot size. Figure 4 shows data obtained in a cross sectional area 
transverse to the direction of ion movement and normal to the surface of the target. The data 
cross section shows ion density versus ion energy level similar to Figure 3; and two 
conditions are shown, generation of the plume under vacuum (4A) and generation of the 
plume in a vacuum but with a partial pressure of nitrogen background gas (Figure 4B). The 
laser pulse was supplied by an Nd:YAG laser system. The experimental system was 
otherwise as per Figure 1. 
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[0028] Figure 5 contains two plots labeled as 5A and 5B. The data of Figure 5 
was obtained by ablation of boron nitride target at 230 picoseconds pulse width, 780 
nanometers wavelength and 6.4 J/cm^ energy fluence. Figure 5 shows ion density similar to 
Figures 3 and 4, and using two conditions, in a vacuum (5A) and in a vacuum with nitrogen 
background gas (SB). The laser pulse was supplied by a TirSapphire laser system. The 
experimental system was otherwise as per Figure 1. 

[0029] Figure 6 shows plots of the ratio of B^"/B^^ for the data as respectively 
shown in Figures 6A (femtosecond case as per 3A), 6B (picosecond case as per 5A) and 6C 
(nanosecond case as per 4A). 

[0030] Figure 7 is a graph showing the signal from the microchannel plate 
(MCP) of the ion energy analyzer in units of voltage as a function of arrival time after the 
laser pulse in microseconds. Figure 7 shows individual charge states and kinetic energies 
for the ionic isotopes of gallium obtained from a gallium nitride target in a vacuum. The 
data of Figure 7 was captured in the experimental apparatus as shown in Figure 1. The 
laser pulse was supplied by a Ti: Sapphire CPA system. 

[0031] Figure 8 shows the charge state and energy distribution data for gallium 
ions and its isotopes, as compiled from spectra such as that on Figure 7. Figure 8 has two 
parts. Figure 8A shows data obtained in the experimental apparatus of Figure 1 where the 
deposition chamber was under vacuum, and Figure 8B shows data obtained in the same 
experimental configuration and under the same conditions except that nitrogen gas was 
included in the chamber. In Figure 8 for the ion species of gallium 69 and gallium 71, there 
is shown data obtained from a cross sectional area transverse to the direction of ion 
movement and normal to the surface of the target. The data shows ion density in ions per 
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square centimeter as a function of kinetic energy expressed in keV. The data of Figures 7 
and 8 were obtained under the conditions of gallium nitride ablation with laser pulses of 150 
femtoseconds, 780 nanometers, 3.3 joules per square centimeter. 

[0032] Figure 9 shows plots of the ratios of Ga^^/Ga''' for three cases. Figures 
9A and 9B are at 150 femtoseconds, 3.3 J/cm^ in vacuum and vacuum plus background 
nitrogen, respectively. Figure 9C is for vacuum at 230 picoseconds, 2.4 J/cm^. 

[0033] Figure 10 is a schematic of a system for depositing a thin film from the 
ablation plume, containing activated, neutral, and ionized particles, onto a substrate. The 
configuration of Figure 10 is similar to that of Figure 1 except in Figure 10 the aperture to 
the detector is blocked by the deposition substrate. As shown in the schematic, a laser beam 
from an external laser source enters the deposition chamber and impinges on a target, 
generates a plume, and the plume is deposited on a substrate holder. 

[0034] Figure 11 shows a schematic illustration of a plume being emitted from a 
target and captured on a substrate. The figure shows that the plume has a Cos°(9) spatial 
distribution with the largest concentration directed normal to the target surface and a 
decreasing concentration away from normal. In the expression Cos"(9), 9 is zero along the 
normal direction. The n represents a power function factor and theta represents the angle 
measured from the direction normal to the surface. With such a distribution, the plume 
deposits a graded concentration of material consisting of a thick portion of deposition on the 
central part of the substrate and a relatively thinner zone of deposition on outer (high 0) 
regions of the substrate. Figure 11 is a schematic representation of the data showing that the 
ratio of B(10) to B(ll) in a central portion of the plume is greater than natural abundance. 
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[0035] Figure 12 shows the Rutherford backscattering spectra of a deposition 
process conducted in accordance with the experimental system shown in Figure 10. Figure 

12 contains results for two samples, one designated as curve 6 and the other as curve 7. 
Curve 6 is for a relatively thin deposition zone of the sample and curve 7 is for a sample 
taken from a relatively thick zone of the deposited film (Figure 11). Both sample spots of 
the film are on the order of 3 millimeters in diameter. The actual spots being analyzed on 
the film were separated by 5 spot diameters or 15 mm (i.e., 1.5 cm). 

[0036] Figure 13 shows the Rutherford backscattering spectrum of a deposition 
process conducted in accordance with the experimental system shown in Figure 10. Figure 

13 contains data for a sample taken from a relatively thick zone of the deposited film. The 
sample spot on the film was on the order of 3 millimeters in diameter. 

[0037] Figure 14 consists of parts A, B and C. Figure 14A shows 5 sampling 
spots of the film taken at separate areas over a 1" region of the film. Sample 1 is at the 
thickest region (plume center) and samples 2-5 are at progressively thinner regions away 
firom the center. The ratio of B^°/B^^ decreases from sample 1 (at B^°/B^^ of 0.45) to sample 
5 (at B^W^ of 0.29). Sample 5 (thinnest region) has B^W^ ratio closest to that which 
naturally occurs (0.25). Figure 14 charts B, C are RBS spectra similar to that shown in 
Figures 12 and 13 but with the Si background subtracted out. The solid lines are 
mathematical best fit curves based on the conmiercially available RUMP computer program 
available through Cornell University. Figure 14B is for spot 5 (thin) and Figure 14C is for 
spot 1 (thick). Note that the vertical scales in 14B and 14C are different. 

[0038] Figure 15 shows isotope ratios, normalized to the natural abundance ratio 
(shown in parenthesis) for boron (15a), zinc (15b), titanium (15c), gallium (15d), and copper 
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(15e) as a function of charge state. Ablation plumes were generated from solid targets with 
780 nm, 200 fs, 2 x 10^"^ W/cm^ laser pulses, to produce results as shown in Figure 15, parts 
(a) to (e). 

[0039] Figure 16 shows angular dependence of boron 1 1 charge state +3 energy 
distribution for ablation of a solid target with 780 nm, 60 fs, 5 x lO*"* W/cm^ laser pulses. 
Angles are measured with respect to the target normal. 

[0040] Figure 17 shows isotopic ratios for boron as observed on the normal axis 
in various laser plasma plumes. Ratios are for BlO/Bll. Natural abundance ratio (0.25) is 
normalized to unity. Pulse durations, wavelengths and intensities are: (a) at 200 fs (NIR), 1 
X 10^* wW; (b) at 6 ns (JR.), 1 x 10^^ W/cm^; (c) at 230 ps (NIR), 3 x 10^° wW; and (d) 
at 6 ns (UV), 4 x 10^ W/cm^; and where NIR = 780 nm, IR = 1.06 /Am, UV = 353 nm. 

[0041] Figure 18 shows a multi-stage arrangement of laser beams and 
enrichment targets. The system is useable with laser beams operated simultaneously or 
sequentially. 

[0042] Figure 19 contains two graphs showing ionic (ni) and neutral (nj 
concentrations in a laser ablation plume formed by 100 fs pulses in Figure 19A and 10ns 
pulses in Figure 19B, at intensities shown in the data box. 

[0043] Figure 20 is a flow chart illustrating one embodiment of the present 
invention in which an initial laser pulse is directed at a target to create a plasma and then a 
time delayed second laser pulse is directed at the plasma to increase the ionic component. 

[0044] Figure 21 is graph illustrating average ion yield and energy as a function 
of time-delay between two identical 120 femtosecond ablation pulses on silicon. The single 
pulse at zero delay has an energy fluence of 2.2 kJ/cm^ on a beam spot diameter of 42 

12 



Attorney Docket No. 2115D-002245 

microns. The two double pulses have a fluence of 1.1 kJ/cm^ each. Expansion distance 
based on measured average ion velocity of 1.9xl0^cm/s 

[0045] Figure 22 is a graph illustrating enhancement of isotope separation for 
boron ions in an ultrafast laser ablation plume. Single pulse: 2.2 kJ/cm"^. Double pulse: 1.1 
kJ/cm'^ each pulse, separated by 10 ps. Laser pulses are 120 fs, 780 nm at 10 Hz. Total 
laser intensity: 2x10^^ W/cm^. Natural abundance = 1. 

[0046] Figure 23 is a flow chart illustrating a second embodiment of the present 
invention in which an initial laser pulse is directed at a target to create a plasma and then 
several individually time delayed subsequent laser pulses are directed at the plasma to 
increase the ionic component. 

[0047] Figure 24 is a graph illustrating elemental enrichment of Ni and Cu as a 
function of angle in the ablation plume. The total enrichment ratio is shown first for each 
angle, followed by each of the ten charge states. The circles and solid line indicate the total 
enrichment ratio. The actual yield of each charge state and angle is not shown in the graph. 
It should be noted that the sample ratio of 45% Ni to 55% Cu was normalized to 1. 

[0048] Figure 25 is a graph illustrating elemental enrichment of strontium (Sr) 
and titanium (Ti) as a function of angle in the ablation plume. 

[0049] Figure 26 includes four graphs showing ion yield as a function of time 
delay between two identical pulses under different laser pulse intensities represented by a 
percentage of maximum intensity (IxlO^^W/cm^) shown on sub-figures a-d. Data points at 
zero delay are produced by a single pulse with twice the energy as each individual time- 
delayed pulse. Maximum laser intensity results are shown as a dash line on the subsequent 
laser intensity results. 
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[0050] Figure 27 includes two graphs showing average ion yield and energy as a 
function of delay time. The intensity of ablation laser pulse is 6% of the maximum 
(6xlO^*W/cm^). Two-component energy distribution is shown in the inset and ion yield and 
energy of each component are also plotted. 

[0051] Figure 28 is a graph showing the charge-state distribution in a silicon 
ablation plume as a function of single full-energy pulse and time delayed (5ps and lOps 
respectively) half-energy pulses. Ion energy distribution is shown in the inset and the peak is 
at 9keV. 

[0052] Figure 29 illustrates a coordinate transformation showing the relationship 
between the one-dimensional cylindrical zone of the plasma centrifuge and the three- 
dimensional expansion zone with the angular distribution associated with the ablation 
plasma. 

[0053] Figure 30 shows an ion energy density distribution for ^^Ni ions as a 
function of charge state. 

[0054] Figure 31 illustrates the angular distribution for the isotope ratio of 58Ni 
to 60Ni as a function of ionic charge state. 

[0055] Figure 32 illustrates the angular distribution of ionic yield as a function 
of charge state for the combined 58Ni and 60Ni ions. 

Detailed Description of the Preferred Embodiments 

[0056] The invention provides new methods for separating isotopes of an 
element and causing enrichment of a desired isotope of an element utilizing lasers to modify 
or fabricate a material containing such isotopes. This invention may be used for a wide 
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variety of materials which contain elements having different isotopes. Before further 
describing the invention, it is useful to understand the problem to which the invention is 
directed. In nature, elements such as boron, gallium, magnesium, carbon and many others, 
occur in stable isotopic form. For example, boron has two isotopes; boron 10, B(10), and 
boron 11, B(ll). Nuclei with the same number of protons but different number of neutrons 
are called isotopes of each other. Therefore, isotopes of a given element have the same 
nuclear charge (protons), but the differing number of neutrons results in a differing atomic 
mass. For a given element, the amount of an isotope in its naturally occurring mixture of 
isotopes, is represented by a term designated as percent natural abundance. These natural 
abundance values are reported in a variety of Standard Chemistry and Physics Handbooks 
and other sources. Using boron as an example, it is known that boron 10 has a natural 
abundance on the order of 19.78% and boron 11 is about 80.22%. It is often desired to 
obtain a material enriched with one of the isotopes relative to the condition found in the 
natural state. 

[0057] In one aspect, the present invention utilizes laser heating and other forms 
of laser induced breakdown of a material to cause ablation and formation of a plasma by 
means of a pulsed laser beam. Pulsed rather than continuous beams may be more effective 
for many procedures, including the formation of an ablation plume from a material 
containing isotopes desired to be separated. A pulsed laser beam comprises bursts or pulses 
of light which are of very short duration, for example, on the order of nanoseconds in 
duration or less. Tj^ically, these pulses are separated by periods of quiescence. The peak 
power of each pulse is relatively high, capable of intensity on the order of 10^ to 10^^ 
watts/cm^. The present invention utilizes one or more laser pulses having selected pulse 

15 



Attorney Docket No. 2115D-002245 

width, repetition rate, wavelength and intensity of the laser beam at the plasma-forming 
material, to form a plasma of the material. Under certain circumstances, continuous laser 
beams are also thought to be useful for the above purposes. Such continuous wave (CW) 
lasers are preferably operated in a range of 10 microns to 200 nanometer wavelength and 
average energy in the range of joules to hundreds of joules. 

[0058] The material containing the isotopes desired to be separated is placed in 
the beam path and is called a target. A plasma is formed from the target material by focusing 
or otherwise dehvering an intense, optical pulse from the laser onto the target. The 
concentrated energy contained in the focused laser beam is absorbed by the target material, 
energizing its electrons and raising it to a temperature or other form of excited state whereby 
ions are produced thereby generating a plasma. The extent of material affected by the laser 
pulse is dependent upon a number of factors including pulse duration, energy, and 
wavelength. 

[0059] The invention provides a method and system apparatus, including laser, 
for producing ions having varying isotopic, ion energy and charge state distributions. In one 
aspect, the ions are contained in a plasma plume generated when one or more laser pulses are 
directed to a target. In another aspect, the ions are collected on a substrate in a manner 
which provides zones having isotopic distribution different from that which naturally occurs 
for a given element contained in the target. It is preferred that the pulses which generate the 
plume be produced by an ultrafast laser (femtosecond or picosecond) which generates the 
plume by ablation or laser induced breakdown (LIB) of the target material. Such ultrafast 
pulses, however, are not thought to be the only way to optimize the process. 
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[0060] Thus, the process and system of the invention produce an unusual isotopic 
enhancement effect in the observed ion spectra. This effect is manifested as an enrichment 
of the lighter isotope in zones of the ablation plume, as observed normal to the surface of the 
target material. The relative intensity of this effect was observed to vary with the charge- 
state and energy of the ions. 

Generation and Analysis of the Plume 

[0061] The process of the invention has been utilized to produce boron and 
gallium isotopes in respective ablation plumes coming from boron nitride and gallium nitride 
solid target. The plume was generated in a system as schematically shown in the diagram of 
Figure 1. As shown in Figure 1, a laser beam was supplied to the deposition chamber 
which contained the target. In the experiments described herein, the laser beam was 
generated by one of two types of laser systems. One was an ultrafast Ti: Sapphire-based laser 
system developed at the Center for Ultrafast Optical Science (CUOS) for the femtosecond 
and picosecond data sets. The other one was a commercial Nd:YAG Laser for the 
nanosecond data set. 

[0062] Nominally, the ultrafast laser system produces 100 femtosecond amplified 
laser pulses at 10 Hertz and a wavelength of 780 nanometers. The pulse energy can be as 
high as 100 millijoules and with a peak-to-background contrast as high as 10^. The laser 
system consists of a number of components and utilizes chirped-pulse amplification (CPA). 
CPA avoids damage in amplifier media by reducing the peak intensity during amplification. 
The ultrashort pulse is stretched in time by chirping the pulse, that is, giving the pulse a 
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time-varying wavelength. The stretched pulse is amplified and then recompressed to give an 
amplified, ultrashort pulse. 

[0063] The oscillator is a commercial Ti:Sapphire system from Clark-MXR, Inc. 
which is pumped by a CW (continuous-wave) Argon laser, also commercial. The mode- 
locked oscillator produces 50-femtosecond pulses at 100 MegaHertz with a pulse energy of 2 
nanojoules. The bandwidth and center wavelength are tunable and are nominally 15 
nanometers and 780 nanometers, respectively. A Pockels cell is used to reduce the repetition 
rate to 10 Hertz which is the (nominal) maximum rate this particular system can handle due 
to thermal effects. 

[0064] The oscillator pulse is stretched in time to 230 picoseconds by a grating- 
based stretcher. The efficiency of the stretcher is 50% and its output is a 1-nanojoule, 230- 
picosecond pulse. The amplification chain consists of three amplifiers. (1) The regenerative 
amplifier is a TirSapphire laser cavity pumped by 5 millijoules of the doubled frequency of a 
Nd:YAG laser (532 nanometers). It can amplify the pulse nominally from 1 nanojoule to 1 
millijoule. (2) The two-pass amplifier is a Ti:Sapphire crystal pumped by 120 millijoules of 
the same Nd:YAG laser used for the regenerative amplifier. Two passes in this crystal 
amplifies the pulse from 1 millijoule to 20 millijoules. (3) The four-pass amplifier is a 
Ti: Sapphire crystal pumped by two high-power Nd:YAG lasers (at a doubled frequency of 
532 nanometers). The total pump energy is 1.2 joules and can amplify the pulse from 20 
millijoules to 200 millijoules. The pulse is then recompressed by a grating pair to a nominal 
pulse width of 100 femtoseconds. The compressor is in a vacuum chamber to avoid self- 
phase modulation and other nonUnear effects of the high-intensity pulse in air. The 
compressor has an efficiency of 50% and outputs 100-femtosecond, 100-millijoule pulses. 
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Background concerning CPA laser systems and lasing using short optical pulse can be found 
in U.S. Patent Nos. 5,235,606 and 5,656,186, incorporated herein by reference in their 
entirety. 

[0065] In the experiments described here, two outputs of the above-described 
laser system were used: (1) the ultrashort pulse which is 780 nanometers and has a dviration 
of 150 or 200 femtoseconds; and (2) the stretched, amplified pulse which is 780 nanometers 
and has a duration of 230 picoseconds. 

[0066] The nanosecond laser used in this research is a Continuum Surelite I 
Nd:YAG laser which has a nominal pulse width of 6 nanoseconds at a wavelength of 1.06 
microns and a pulse energy of 450 millijoules. 

[0067] For the experiments, the beam had pulse widths in the nanosecond (ns) to 
femtosecond (fs) range, repetition rate on the order of 10 Hz, intensity in the range of 10^° to 
10^^ watts per square centimeter (w/cm^), and wavelength on the order of one micron. The 
laser beam was focused and directed to the target through a lens. 

[0068] The laser pulse was directed at the target at an angle to the target of 
approximately 45 degrees. The laser beam ablated the target material and generated a 
plume. The specific angle of laser impingement on the target is not thought to be critical. 
However, it is important to note that the ablation plume leaves the surface of the material in 
a direction that is essentially normal to the surface, that is, at a 90 degree angle to the 
surface. The target could be any material including amorphous or crystalline solids, liquids 
or gases including high density plasma/gas state. The examples described herein were 
solids. 
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[0069] The center of the plume was directed toward a set of apertures that acted 
as collimating devices and introduced the ions into the ion energy analyzer. This analyzer 
was a conventional curved (spherical) sector plate electrostatic energy analyzer. When it 
was operated in a time of flight mode, through synchronization with the laser pulse, the 
analyzer was also capable of uniquely separating ion charge states. Discrete energy 
distributions of the individual charge states for ions in the ablation plume were obtained by 
sweeping the analyzing voltage on the sector plates and using time of flight synchronization. 

Monitoring Ion Charge, Energy and 
Distribution of Charge and Energy 

[0070] The energies and charge-states in the ablation plume reported herein were 
measured using a spherical-sector electrostatic energy analyzer. The ions entered the 
analyzer and passed between two curved spherical-sector plates. A voltage up to 5000 volts 
was placed on each of the plates. An ion passing between the plates is influenced by this 
voltage which curves the path of the ion. The magnitude of this force depends on the voltage 
on the plates and on the charge-state and energy of the ion. Specifically, for a certain voltage 
difference on the plates, only ions with a certain energy-to-charge ratio will follow the path 
and reach the detector (the MCPs). In the geometry of this analyzer, the energy-to-charge 
ratio (E/q) of the selected ions is equal to 2.254 times the voltage difference (Vdif?) between 
the spherical-sector plates which gives the equation: (E/q)seiected = 2.254 x Vdiff where E/q is 
in units of electron volts (eV), and Vdiff is in volts (V). For positive ions, a negative voltage 
was placed on the inner plate and a positive voltage on the outer plate. 
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[0071] A specific voltage difference between the plates gives ions with a 
particular energy-to-charge ratio. The detector will observe a charge state 1 ion with a 
specific energy E as determined by the sector plate voltage and will also observe a charge 
state 2 ion at an energy of 2E. However, ions with differing charge-states are detected at 
differing times after the laser pulse. This follows since an ion with twice the charge-state 
and twice the energy travels at a faster speed, that is, 1.414 times faster because energy is 
related to the square of the speed. Thus, the higher-energy ion arrives earlier. This 
technique which used the arrival time after the laser pulse, to analyze the detected signal, is 
called spectral time-of -flight analysis. 

[0072] In order to obtain the entire distribution of ions over energy and 
charge-state, the voltage difference on the sector plates was scanned from 0 to 10,000 volts. 
Each voltage difference gave a set of charge-states and energies as described above. A 
computer program is used to select the individual charge-states at each sector plate voltage 
difference and to assemble the ion energy and charge-state distributions. These distributions 
of ion energies and charge-states depend on the parameters of the laser pulses as defined 
herein. The isotope enrichment process was observed both in the ion distributions of the 
plume and in the spatial patterns, from the deposition and growth process, on the thin films. 

[0073] A typical spectrum is shown in Figure 2 for boron nitride where the 
boron isotopes are shown along with the nitrogen ions. For the case of a boron nitride target, 
the two stable isotopes, boron 10 and boron 11, were easily observed. Their natural 
abundance is 19.78% and 80.22%, respectively. Relative isotopic ratios in the ablated plume 
are observable in these spectra as a function of ion charge state and ion energy. 
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[0074] Figure 2 shows data obtained at a fixed voltage on the sector plates and 

using time of flight synchronization to build up discrete energy spectra of the individual 
charge states for ions in the ablation plume. Here, the signal from the MCP (micro channel 
plates) of the detector was plotted as the function of arrival time after the laser pulse shown 
in microseconds. This data was obtained under vacuum conditions. This data was obtained 
using pulse width of 200 femtoseconds, wavelength of 780 nanometers, repetition rate of 10 
Hz, and 4 millijoules (mJ) energy per pulse using the Ti:Sapphire laser. The data of Figure 
2 was captured in a vacuum by the experimental apparatus shown in Figure 1. The intensity 
of the laser beam incident on the target was about 2.5 x lO'"* watts per square centimeter. 

[0075] Referring to Figure 2 and bearing in mind that the natural abundance 
ratio of boron 10 to 11 is on the order of 20:80, it can be seen that enrichment of a specific 
isotopic species within a zone of the plume is achieved. More specifically. Figure 2 shows 
boron 10 B(10) and boron 11 B(ll) at the +4 charge condition captured approximately 4 
microseconds after the laser pulse. Here, the occurrence of B(10) and B(ll) is about the 
same. In other words, at the +4 charge condition corresponding to 4 keV, the number of 
B(10) isotopes, as shown by the B(10) peak, is approximately the same as the number of 
B(ll) isotopes. This distribution of B(10) and B(ll) at the +4 charge state shows an 
occurrence of B(10) that is 4 times more than its naturally occurring abundance with respect 
to B(l 1). Therefore, there was clearly a zone within the plume that contained an enriched 
fraction of B(10) with respect to B(l 1) as compared to the naturally occurring amount. The 
distribution captured for B(10) and B(ll) in the +1 charge state at about 1 keV, is shown in 
the vicinity of 8-9 microseconds after arrival of the laser pulse. In the graph at close to 9 
microseconds, there is shown a peak for B(l 1) and close to 8 microseconds, there is shown a 

22 



Attorney Docket No. 2115D-002245 
peak for B(10). The relative proportion of B(10) to B(l 1) at this segment shows that there is 
more B(10) relative to B(ll) than the naturally occurring amount. This is even more clearly 
seen in Figure 3A, where the ratio of B(11):B(10) is always less than 4:1 and often at about 
1.5:1. The naturally occurring amount of B(ll) to B(10) is 4 to 1 (80:20). This enrichment 
at essentially all charge states +1 to +5 occurred because the aperture permitted passage of 
the rich B(10) part of the plume to the detector. Correspondingly, B(10) lean parts of the 
plume did not reach the detector. Since isotopic species were not created or destroyed in the 
process, it follows that there are probably zones within the plume or remaining target 
material which are relatively deficient in B(10) and relatively enriched in B(l 1). 

[0076] The data of Figure 3 was obtained by sweeping the voltage on the 
analyzed sector plates, the charge and energy distribution data for boron as shown in Figure 
2 was integrated to produce Figure 3A. Figure 3 also shows data for equivalent 
experimental conditions where a deposition chamber included nitrogen gas backfill of 0.8 m 
Torr (Figure 3B). There are 3 sets of data in Figures 3A and 3B, respectively, showing 
B(ll), B(10) and N (nitrogen). The nitrogen curves are significantly smaller in amplitude 
than the boron curves due to the depletion of N in the BN target from heating and ablation 
over time. 

[0077] Figure 3A shows the ion energy distributions resulting from successive 
spectra as obtained per Figure 2. In Figure 3A the peak of the +1 ion energy distribution is 
to the extreme left followed by +2, +3, +4, and +5. B(ll) is the top plot and B(10) the 
middle plot. Peaks from nitrogen are also given as the bottom plot of Figure 3A, but are not 
being used to demonstrate the effect in this figure. It is anticipated that the very small 
concentrations of nitrogen isotopes, not desirable here, are likewise undergoing enrichment. 
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These nitrogen peaks exist because the starting material is boron nitride and nitrogen ions 
were generated as well as boron ions. More specifically, Figure 3A shows a plot of 
transverse ion density expressed in ions/cm^ versus ion energy expressed in keV. Figure 3A 
was obtained by sweeping the sector voltages on the analyzer plates resulting in the energy 
distributions as a function of the charge state. As can be seen from Figure 3A, the geometry 
of the experiment of Figure 1, caused capture and detection of the zone within the laser 
ablation plumes where the relative ratio of these isotopes (B(10) and B(ll)) are drastically 
altered from their naturally occurring abundance, when measured as a function of charge 
state and energy. That is, in all cases, the ratio of B(10):B(1 1) was greater than that which 
naturally occurs, and often close to 50:50. 

[0078] The results demonstrated in Figures 2 and 3A were also obtained when 
generating a plume in the deposition chamber in the presence of a background gas of 
nitrogen. When the chamber contained nitrogen gas, there was a charge exchange which 
occurred as the boron ions gathered electrons from the nitrogen. Here, the boron ions were 
reduced from the higher charge states such as +3, +4 and +5 and to the lower charge states, 
namely, +1 and +2. This is evident from Figure 3B. Here, enhancement was also clearly 
obtainable since the ratio of B(10) to B(l 1) was on the order of 50:50. 

[0079] It is thought that the presence of the background gas facilitates formation 
of a product film containing constituents derived from the background gas and this may be 
desirable in forming isotope enriched zones having a particular chemical composition. It is 
clearly advantageous that the presence of a background constituent does not affect the 
process of enrichment yet permits inclusion of background gas species in a product. 
Therefore, by this means, materials of desired composition may be formed having the 
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isotope enriched feature. This is demonstrated later with respect to formation of a film from 
the plume. It is anticipated that isotopically enriched particulates, precipitated from the 
plume or background gas, could also be formed. This would provide a mechanism for 
creating isotopically enriched nano-particles. 

[0080] Figure 4 contains results obtained from the experimental apparatus of 
Figure 1 but with a different laser beam. Here, laser wavelength of 1.06 microns was used; 
repetition rate of 10 Hz was used; and a pulse width of 6 nanoseconds was used. The pulse 
energy was 60 millijoules (mJ) energy per pulse corresponding to a fluence of 760 J/cm^ and 
an intensity of 1.3 x 10^^ w/cm^. Here, the main difference in the method was the use of 
significantly different pulse width, namely, 6 nanoseconds as compared to the prior 200 
femtoseconds. The experiment was conducted both in a vacuum and in the presence of 
nitrogen background gas. 

[0081] The experimental set up used to obtain data of Figure 4 included the 
directing of the plume through an aperture as shown in Figure 1. Therefore, the data of 
Figure 4 also shows capture and detection of the B(10) enriched central portion (cone) of the 
plume. 

[0082] Figure 4 shows transverse ion density of the plume captured in the 
enriched zone as a function of energy level. The results shown in Figure 4 are both for the 
case of the vacuum and the case where nitrogen background gas is contained in the 
deposition chamber. It is observed that generally the transverse boron ion density in the 6 

nanosecond case is less than the 200 femtosecond case. In other words, there is less dense 
plasma plume generated under the laser operating parameters defined here at 6 nanoseconds 
pulse width. Relative enrichment is still obtainable. In other words, there is enrichment of 
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B(10) with respect to B(ll) as compared to that which is naturally occurring. The 
overwhelming conclusion to be drawn from the data is that enrichment has clearly occurred. 
The effect of N2 background gas in Figure 4B was similar to that observed in Figure 3B. 
Therefore, the overall conclusion to be drawn from the data presented in Figures 2, 3, and 4 
is that the relative ratio of the isotopes, in these ablation plumes, are drastically different 
from that observed for the naturally occurring species as reported in standard references. 
More specifically, the experimental set-up shown in Figure 1 contains the apertures having 
four millimeter (mm) holes. The aperture guided the plume into the detector and also acted 
as a way of isolating the high pressure zone of the chamber from the low pressure (vacuum) 
zone. Since the plume was generated in a vacuum chamber, apertures help in maintaining 
the relative vacuum. As a consequence of this set-up, the aperture essentially permitted 
passage through of that portion of the plume which contained the relatively high proportion 
of B(10) as compared to B(ll). Therefore, all of the data shown in Figures 2, 3, and 4 show 
a relatively enriched B(10) portion of the plume. 

[0083] Figures 5A and 5B show the isotope enhancement effect for a picosecond 
laser pulse. The energy and charge state distribution of the ions are shown from a BN 
plume. The plume was generated with a 780 nanometer, 230 picosecond laser pulse with an 
intensity of 2.8 x 10^° Watts/cm^. Figure 5A is for a plume in vacuum, and Figure 5B is for 
the same plume in a background of 0.8 milliTorr of N2 gas. Charge state +3 in the vacuum 
case and charge state +1 in the N2 case show a very strong enhancement of B(10) over 
B(ll). This demonstrates the dependence of the isotope effect on ion charge state and ion 
energy. 
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[0084] The dependence of the isotope enhancement on ion charge state, ion 
energy and laser pulse duration is summarized in Figure 6. The calculation of the B(10) to 
B(ll) ratio for the data from Figure 3A, 5A, and 4A is shown in Figures 6A-C, 
respectively. The naturally occurring ratio is 0.25 and is shown as the dotted line in each 
graph. The femtosecond case in Figure 6A shows a strong enhancement of B(10) over 
B(l 1) for charge states +1 to +5, particularly for +3, +4, and +5. The enhancement is present 
for nearly the entire range of the ion energies. 

[0085] The isotope enhancement was also measured for gallium which has 
naturally-occurring isotopes with masses of 69 and 71 amu with an abundance of 60.4% 
Ga69 and 39.6% Ga71. Enrichment of the lighter isotope was also observed for this material 
using the experimental set-up of Figure 1 in the same manner as for boron nitride. The 
time-of-flight spectrum from a GaN (gallium nitride) ablation plume in vacuum is shown in 
Figure 7. The laser pulse was 150 femtoseconds, 780 nanometers, and 2.2 x 10^^ 
Watts/cm^. In charge states +1 to +3, Ga69 and Ga71 are present, with ratios clearly greater 
than the naturally-occurring ratio of 1.53. 

[0086] The time-of-flight spectra are compiled in Figure 8 for the same laser 
pulse described for Figure 7. Figure 8A shows the ion charge state and energy distribution 
in vacuum and Figure SB shows the case in 0.8 milliTorr N2 background gas. Once again, 
the enhancement of the lighter isotope, Ga69, over the heavier, Ga71, above the naturally- 
occurring abundance is clear for the different charge states and energies. 

[0087] The isotope ratio for gallium is shown in Figure 9 for three conditions. 
Figures 9A and 9B are for 150 femtosecond, 780 nanometer, 3.3 Joules/cm^ laser ablation 
plumes. Figure 9A shows the calculated Ga69 to Ga71 ratio for the plume in vacuum, and 
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Figure 9B in a background gas of 0.8 milliTorr N2. Figure 9C shows the ratio for a 230 
picosecond, 780 nanometer, 2.4 Joules/cm^ laser ablation plume in vacuum. The naturally 
occurring ratio 1.53 is shown as dotted lines in the figures. The enhancement of the hghter 
isotope, Ga69, over the heavier one, Ga71, follows the same trend as the boron case for 
B(10) and B(l 1). The effect is generally stronger here for gallium than that shown for boron 
in Figures 3-5. 

[0088] The effect observed in the time of flight ablation plume coming from the 
ablated target is useful for preparing film products. This is accomplished by directing the 
plume toward a substrate in an experimental arrangement as shown in Figure 10. As shown 
schematically in Figure 11, the plume condenses to form a thin film on the surface of the 
substrate. The film was enriched in the preferential isotope ratios as observed earlier and as 
demonstrated by Figures 2-6. The following film experiments were designed to capture 
essentially the entire cross section of the plume. By this means, zones rich and lean in B(10) 
as compared to the natural abundance were clearly observed and verified. More specifically, 
it is believed that isotopes are neither created nor destroyed in the process, therefore, it is 
assumed that the composition of the plume or remaining affected target material, in total, 
represent natural occurrence of the B(10) and B(ll) isotopes. This is further illustrated by 
reference to Figure 11. Figure 11 is a schematic showing the plume coming off of the 
target. The plume has a relatively B(10) enriched conical section (zone) and a relatively 
lesser B(10) zone on either side of the concentrated zone. These lesser concentration zones 
may be B(10) depleted. In order to further understand and verify the relative proportions of 
the isotopes in the entirety of the plume, films were formed in the deposition chamber in a 
manner which captured essentially the entire cross-section of the plimie. This entire cross- 
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section encompassed a relatively B(10)-rich central conical section as well as surrounding 
zones having relatively less B(10) or being B(10) depleted. 

Film Production and Analysis 

[0089] Figure 10 shows an experimental set-up very much similar to the 
experimental set-up of Figure 1 except that the plume generated from the target is directed 
to a substrate to form a film thereon. No apertures are used to screen out any portion of the 
plume or film. As shown in Figure 10, the laser generates a pulse of 150-200 femtoseconds, 
780 nanometers, and 10 Hz with an energy per pulse of 4mJ-100mJ. The laser beams are 
directed through a focusing device onto the target. In this example, the target is boron 
nitride. A plume of laser plasma is generated containing ionized particles of boron, nitrogen, 
and also electrons, some neutrals, and likely some excited atoms and molecules. The plume 
is directed onto a 2-inch diameter silicon substrate where the film is formed. Optionally, the 
substrate holder may be cooled by hquid nitrogen to a temperature on the order of 77°K if 
desired to affect the condensation and morphology of the formed film. Substrate 
temperatures elevated above room temperature are also possible up to 1000°C. The 
electrostatic analyzer (charge state and energy analyzer) is also shown. 

[0090] The silicon substrate had a thickness on the order of 400 microns and 
boron nitride films were deposited thereon having a thickness of approximately 0.1 microns. 
Once the system had been set-up, aligned and calibrated, films were produced from the 
plasma in a relatively short time frame on the order of about 15 minutes. 

[0091] The results of the deposition process are shown in Table I and 
corresponding graphs of Figures 12 and 13. Table I contains as a first entry, a listing of the 
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relative abundance in nature of B(l 1) and B(10). The natural occurring ratio of B(l 1):B(10) 
is about 4:1. The data entry designated as BN3 169906 (-906) was taken from a relatively 
thin zone of the film. As shown in Figure 11, the film thickness decreases from the center 
toward the outer edge of the substrate. An analysis of this thin part of the film, based on 
Rutherford ion beam analysis, shows the counts for B(10) and B(ll). The ratio of B(I1) to 
B(10) at the -906 data point is relatively rich in B(ll) (i.e., depleted in B(10)). Here, the 
amount of B(10) is on the order of 17% whereas its naturally occurring amount is close to 
20%. The amount of B(ll) in the -906 sample is close to 82%, whereas its naturally 
occurring amount is about 80%. 

[0092] Samples 3169907 (-907) and 3169908 (-908) were taken from portions of 
the film designated as thick and near the center of the substrate. This is also schematically 
shown in Figure 11. In the -907 sample, B(10) was on the order of 30% and B(ll) on the 
order of 70%. Here it is clearly seen that the amount of B(10) isotope present in this zone of 
the film is very much greater than the naturally occurring amount which is about 20%. 
Sample -908 (Figure 13) confirms the results found in sample -907. In sample -908, isotope 
B(10) was present in an amount close to 35% based on the sample count of B(10) by the 
Rutherford scattering method, and B(l 1) count was on the order of 65%. Here it can be seen 
that the amount of B(10) is very highly enriched as compared to its naturally occurring 
amount with respect to B(ll). Referring back to curve 6 of Figure 12, (sample -906 taken 
from the thin portion of the film) it is possible to subtract out the effect of the silicon 
substrate by virtue of the line, and look above the line for the peak in the vicinity of channel 
200 which shows B(10) in a relatively small proportion for curve 6. The proportion of B(10) 
in the thick portion of the film is the upper curve labeled curve 7. Here the proportion of 
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B(10) is greater than that which naturally occurs. B(l 1) peak is shown closer to the channel 
300 on the curves. Overall, it can be seen that the proportion of B(ll) to B(10) in curve 6 is 
higher than the proportion of B(ll) to B(10) in curve 7. 

[0093] Figure 13 contains results corresponding to sample -908. Sample -908 
was taken near the relatively thick center of the film. B(10) is shown as a peak in the 
vicinity of channel 200 and B(ll) is shown in the vicinity of channel 300. Here the 
proportion of B(10) is higher than either of the samples -907 or -906. It should be noted that 
the peak near the channel 400 of the X axis of all the curves corresponds to nitrogen. In all 
cases, the baseline represents the silicon substrate. 

[0094] The BN film was analyzed by Rutherford backscattering (RBS) at 5 spots 
of 3 mm size over a range of 1", as shown in Figure 14 A. The RBS spectrum was fit using 
the commercially-available computer program RUMP. The silicon background was 
subtracted out and the B(10) to B(ll) ratio was fit as shown in Figures 14B and 14C for 
different spots on the film. The thickest region of the film was sampled at spot (1) which 
corresponds to the deposition from the central part of the plume. The thinnest region of the 
film was sampled at spot (5) which corresponds to the outer part of the plume. The 
B(10)/B(ll) ratio is maximum for the film captured from the plume center and decreases 
toward the naturally-occurring ratio away from the center. 

[0095] The experimental set-up of Figure 1 was used in ftjrther experiments, 
using the method applied to boron and gallium. 

[0096] The method applied to boron and gallium produced measured factors 
greater than 2 for isotopic enrichment in the multiple-charge-state ions from ultrafast laser 
ablation plasmas. Such methods applied to the elements Zn, Ti, and Cu produced results 
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which confirm the enrichment of the lighter isotope on the target-normal axis of the 
expanding ablation plume (see Figure 15). Angular scans of these plumes for boron 
demonstrated the persistence of the enrichment effect to wide angle, being present as far out 
as 45 degrees from the normal direction (central axis). Systematic variations in the dominant 
ion species and their energy distributions, as a function of charge state, are observed in these 
angular scans (see for example Figure 16). 

[0097] Deposition of isotopically-enriched thin films appear to be strongly 
influenced by the ratio of ionic to neutral-atom species in the plasma. The data here 
demonstrate a relationship between the ion to neutral atom ratios in the plumes and the laser 
operating parameters. A variety of techniques were employed to show such ratios as a 
function of laser operating parameters, including deflected versus non-deflected neutral 
species in the ion energy analyzer, laser probe deflection in the refractive portion of the 
ablation plasma, emission spectroscopy, and Langmuir probe measurements. In general, the 
neutral species were observed to have several factors of 10 lower energies than the ions. 
Likewise, the data indicates that laser operating parameters are used and adjusted to 
maximize or minimize the ionic to neutral specie ratios in the plume. It is thought that 
secondary laser pulses which act as plasma pumping pulses could significantly change the 
ion to neutron ratio and thereby provide greater enrichment. 

[0098] Figure 15 shows the isotopic ratios of several elements from the periodic 
table, as a function of charge state, as observed in ultrafast laser ablation plasmas. In Figure 
16 are the results of an angular scan showing the variation in energy distribution, as a 
function of angle from normal for moderate-intensity ablation (5 x lO'"* W/cm^). These 
results are consistent with the model of a plasma centrifuge for separation of ions according 
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to kinetic energy. This plasma centrifuge is thought to be controlled by the longitudinal 
magnetic field of the laser ablation plasma. It is thought that the expanding plasma produces 
a strong magnetic field that provides a rotational component to the ions as they explode 
(expand) away from the ablating surface. 

[0099] These results demonstrate a kinetic energy separation for the ion energy 
distributions as a function of observation angle relative to the target normal. Such kinetic 
energies are proportional to the first power of the mass and second power in velocity (K.E. = 
V2 mv^), thus confirming the presence of lower mass ions (with lower velocity) near the zero 
(or central) axis. 

[0100] The results shown in Figures 15 and 16 for the various elements are 
consistent with the results obtained for boron. Returning to the data for boron, Table II 
shows the integrated areas under each of the charge- state curves presented in Figure 3. The 
ratios of these areas are also presented which provides a measure of the total isotope 
enrichment associated with each charge state of boron. An interesting pattern is seen in the 
case of vacuum ablation where charge state +3 has the largest (nearly 1:1) ratio of the light to 
heavy isotope. This is followed by charge state +4 which has the next most highly-enriched 
condition. It is also observed that the highest total yield of ions is for the +4 charge state. 

[0101] In summary, the invention provides an isotopic enrichment process in 
laser ablation plumes. This effect was demonstrated in the ion energy spectra of laser 
ablated boron nitride and gallium nitride materials. The effect was first observed using 150- 
200 femtosecond pulses and has also been seen using 200 picosecond and 6 nanosecond 
pulses. It appears to be more pronounced in the short pulse case (fs) but is also clearly 
observable in the longer pulse case (ns). The manifestation of this effect occurs in the yields 
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obtained for the boron 10 compared to boron 1 1, as well as for gallium 69 versus gallium 71. 
These isotopes are present in the ionically charged species of the laser ablation plume. The 
naturally occurring abundance of boron 10 to boron 11 is 19.78% to 80.22%, nominally 
20/80 ratio. In the ion energy spectra reported for the enriched zone of the plume, as a 
function of ion-charge state and ion energy, there are ratios that are greater than 50/50 and in 
most circumstances always greater than the naturally occurring 20/80 ratio. These 
differences seem to vary with ion charge state, ion energy and laser pulse duration. In 
addition, the isotope enrichment is generally larger for gallium, which suggests an important 
dependence on the target mass. 

[0102] In the systems described herein, the laser pulse was focused to the surface 
of the ablation target and the ablation plume left the surface of the target approximately in a 
normal direction expanding into the region surrounding it. This region constituted a high 
vacuum or an enclosure containing a dilute background gas. The particles in the plume 
consisted of electrons, ions, and neutral atoms or molecules. These neutrals may or may not 
be in excited states. The plume was captured so as to provide zones containing isotope 
distribution different from the natural abundance. The results, as presented above, are 
unusual and unexpected. This is because isotope separation, as a process, is normally done 
by taking advantage of known features such as mass diffusion or difference in atomic 
vibrational modes. The process described herein is completely different. 

[0103] While not wishing to be held to any particular theory, one possible 
explanation for the effect could be found in the nature of atomic collisions for large 
backward angles in the present case. The laser ablation process is produced when the laser 
beam is absorbed by the material at its surface resulting in an intensely ionized and very 

34 



Attorney Docket No. 2115D-002245 

dense plasma. This plasma expands outward away from the surface and inward at the same 
time. The inward going portion is reflected by the solid surface beneath it and is scattered 
back in a rearward direction reinforcing the original outward moving plasma front. For 
ultrafast pulses (fs and ps) the plasma at the surface will be very thin (a few hundred 
angstroms) and the collisions producing the backward directed motion will consist of 
predominantly ions that scatter at approximately 180 degrees. The observation angle for the 
detector is at exactly 180 degrees with, in addition, a very small acceptance angle, on the 
order of 0.1 degree.In reviewing the above results, several possible phenomena and 
properties are to be considered. They are summarized here. The effect is restricted to the 
first few hundred angstroms below the surface of the material. Screened coulomb collision 
mechanisms and other possible screened potentials seem to pertain. Mass differences play 
an important role in the process. Shallow angle collisions just prior to the main ultra-low 
impact parameter backscattering event seem to be relevant to the process. Recoil of the 
target nucleus seem to be a factor in the enhancement. The angular width of the effect varies 
with ion energy, charge state and target mass. The effective depth over which the process 
occurs is a function of ion mass and energy and the target mass. In addition, variations in the 
plasma distribution can produce a type of electromagnetic lens which selects and focuses 
ions according to charge state and energy. 

[0104] Therefore, a filtering effect can be based upon coUisional filtering or 
electromagnetic filtering. It is clear that the plasma plume has an angular distribution pattern 
with the central cone preferential with the lighter isotope. This is thought to be because 
when the plume comes off the target, there is an angular distribution pattern for the isotopic 
species. 
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[0105] All of these parameters seem to be important when comparing light and 

heavy isotopes that are involved in a wide angle 180 degree scattering event. By using 
scaling, it is expected that the angular width of the effect observed here will be on the order 
of 8 to 10 degrees from the surface normal. This will result in an isotopically enriched 
circular zone of about 1 centimeter radius near the center of the deposited film, if that film is 
deposited at a typical 7.5 centimeters away from the ablation source. At one meter from the 
ablation source the angular width will produce a circular zone of enrichment that is on the 
order of 14 centimeters in radius. It is evident that this process is adaptable to commercial 
scale-up. 

[01 06] As stated above, data was collected for the 6 nanosecond pulse duration, 
150 to 200 femtosecond pulse duration, and also for an intermediate pulse duration of 230 
picoseconds. With gallium, the effect at 230 picoseconds appears to be diminished relative 
to the 150 femtosecond case. With boron, there is no clear trend in B(10)/B(ll) ratios from 
200 femtoseconds to 230 picoseconds to 6 nanoseconds. In particular, there is a charge state 
dependent isotope enhancement for certain conditions, shown most strongly in the +3 charge 
state of Figure 6B. However, the total yield is decreased in the nanosecond case, relative to 
the picosecond and femtosecond cases. The gallium and boron data suggests that there is a 
relative optimum pulse duration for obtaining the ablation effect. 

[0107] As a result there is a very thin layer on the order of 200 angstroms where 
the process occurs, the correlated collisions occur within a very shallow depth. This 
enhancement or optimization of the backscattering mode seems to be dependent on the mass, 
energy and charge. It appears that the cone in which the effect occurs is dependent on the 
mass of the scattering atoms. The heavier the mass of the scattering atoms, the bigger the 
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cone. Conversely, the lighter the mass of the scattering atoms, the smaller the cone. From 
the data obtained in the B(10) and B(ll) relationship, it appears to show that the mass effect 
is there and the charge effect is there. The screen charge appears to be related to why the 
different charge states provide different isotopic ratios. Therefore, the data clearly suggests 
there is a mass dependence in the scattering and a charge dependence in the scattering and 
the width of the cone is controlled by that. Both of these dependencies namely charge and 
mass are observed in this work. 

[0108] Backscattering as used herein designates travehng of particles in a 
direction opposite to that of the path of the energy. In this case, backscattering is not related 
to the angle of incidence of the laser. In contrast, backscattering as defined herein is the 
reaction of the particles to the movement of the plasma in a direction into the material. This 
effect is the result of the formation of the plasma after the laser pulse has delivered its 
energy. 

[0109] It is thought that pulse width on the order of hundreds of picoseconds or 
less, which cause ablation in plasma formation while minimizing thermal effects, are 

preferred. This generates a relatively dense plasma from which the ions expand. There may 
be a slight material dependence with respect to the preferred pulse duration or pulse width. 
Methods for generating short pulses such as 200 ps or less or 10 ps or less are described in 
U.S. Patent No. 5,656,186. 

[0110] It is not fully known from the present experiments if the process is 
significantly wavelength dependent. Preliminary experiments suggest the effect is very 
small, or may not be present, for 353 nm (i.e., near UV). The work was done in the near IR 
at about 1 micron. There could be useful results for the near UV wavelength range since it is 
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known that such laser pulses, even when relatively long and in the nanosecond range, are 
absorbed over a very thin layer of only a few hundred angstroms. It is therefore thought that 
the process of the invention is workable in the ultraviolet wavelength range with longer 
(nanosecond) pulses. 

[0111] Therefore, wavelength on the order of 198 nanometers (nm) to 2.5 
microns is thought to be useful along with pulse widths on the order of nanosecond to 
femtosecond with energy fluences in the range of millijoules/cm'^ to hundreds of joules/cm'^ 
or more, and intensities on the order of 10^ to 10*^ watts/cm^. Extreme relativistic electron 
effects become significant above 10^^ watts/cm^ and could provide new processes relating to 
the herein observed isotope enrichments. It is thought that the preferred range is 0.8 micron 
wavelength with 100 femtosecond pulses at lO'^ to 10^'' watts/cm^ laser intensity. 

[0112] In the case of femtosecond short pulses, the absorption depth of the 
energy is relatively shallow. In the case of a relatively long wavelength, there may be a 
trade-off between a relatively deep absorption depth and the length of the pulse. For 
example, a longer pulse at a longer wavelength is thought to affect a greater depth of 
material and the isotope enrichment effect may be less as per the data in the 6 nanosecond 
case. 

[01 1 3] The process does not have to be done in a vacuum, however that is the 
cleanest way to get a film. It could be done in a background gas (as shown earlier) and even 
conceivably in air. Doing this in air (or other gas at atmosphere, or above) could, in fact, be 
an interesting way to make isotopically enriched particulates. The size of these could be 
controlled by the gas pressure. 
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[0114] It is not thought that repetition rate of the laser plays an important role 
except that the growth rate of the film will be affected thereby. The energy level of the laser 
is important and it is desired to ablate material and form the plume without residual thermal 
effects into a substantial bulk of the ablated material. 

[01 1 5] Energetic ions have a very wide range of initial energies and charge states 
in the laser ablation plasma. As such, they represent a very rich source of diagnostic probe 
particles. Appropriate detectors can examine spatial and energy distributions of various 
charge states in order to understand phenomena taking place within the ablation plasma. 
Experimental results herein show how such swift ionic charge states are affected by passage 
out of the ablation plume and what effect occurs as a function of observation angle about the 
normal to the ablating surface. These experiments make it clear that observing energetic 
ions as a function of energy, mass, and spatial (angular) distribution represents a rich source 
of information. Such data provides valuable information on the force fields affecting these 
ions beyond conventional fluid dynamic forces and hydrodynamic processes. Figure 17 
shows an example of such measurement where the effects of laser wavelength and pulse 
duration are examined relative to isotopic mass ratio analysis. 

[0116] Interesting variations of isotopic ratios are seen as a function of laser 
wavelength and particle charge-state. It is thought that these variations are directly related to 
the time development, strength, and orientation of the magnetic fields in the ablation plasma 
through which the swift ions pass. Ions of different charge state travel at different velocities 
relative to the beginning of the ablation pulse. These ions then pass through the plasma 
magnetic fields at different times relative to the build up and decay of the magnetic fields. 
Thus, such ions serve as a probe for determining the growth rate and duration of such fields. 
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It is noteworthy that the isotopic separation of ions seen in Figure 17 is a minimum for 
experiments done in the UV (353 nm) compared to the IR (1.06 micron) and near IR (780 
nm). This implies a wavelength dependence of the magnetic-field-generating mechanisms. 
The angular distribution of ejected ions as a function of charge state, energy, and mass were 
examined. These experiments suggest strongly that there are pronounced kinetic energy 
dependencies to such angular distributions. Again, this is thought to be due to the presence 
of magnetic fields within the plasma through which the ions penetrated. 

[0117] Figure 18 shows a multi-stage arrangement of laser beams and 
enrichment targets. The system is operable with all laser beams and targets run 
simultaneously or sequentially. In one embodiment, all targets are the same material, and as 
time goes on each target is enriched in a selected ionic species. The last target down the 
chain will have the greatest enrichment, due to the multiplication effect that occurs. If it is 
desired to achieve enrichment of the lighter isotope, then the laser is directed to the inner 
conical region of the deposit formed by ablation of the upstream target in the cascade 
arrangement. It if is desired to achieve enrichment of a heavier isotope, then the laser is 
directed to the annular region outside the inner conical region of the deposit formed by 
ablation of the upstream target. 

[01 1 8] Assuming that the lasers and targets are nominally the same and operation 
of the laser is the same, a dynamic equilibrium condition arises. Here, the rate of material 
added, for example, at target 2 is equal to the rate of material removed from target 1. After 
some desired period of operation, the last target in the chain is the most enriched with 
respect to a selected isotope. 
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[01 19] In another embodiment, the irradiation of targets is conducted sequentially 

to provide progressively enriched films in the cascade arrangement of Figure 18. Thus, an 
initial spatial separation of isotopes occurs by irradiating (ablating) target 1 whereby a 
deposit is formed having a relatively enriched region and lean region for a selected isotope. 
Such deposit forms the second target and one of the regions thereof is irradiated to provide 
further isotope separation. 

[0120] Significantly, the isotopic enrichment is clearly observed in the ablation 
ion spectra and in the deposited films. The commercial implications are evident. Products 
having desired isotope distribution are able to be produced from the process. Multiple 
iterations of film formation may be used to increase the isotope concentration. The film 
formed in a prior iteration becomes the subsequent target to produce a desired result. 
Multiple deposition steps are also useable to form multi-layer films. The process and system 
apparatus are adaptable to commercial use and automated production. Therefore, isotope 
enriched products are obtainable by the methods and system of the invention. 

[0121] Although not wishing to be held to a particular theory, the enrichment 
effects observed in this work most probably involve interaction of charged ion species in the 
laser plasma with intense magnetic fields formed by the laser just ahead of the plume 
expansion process. It is less likely that the observed effects are associated with isotope 
production through laser-induced nuclear reaction processes. The laser intensities used in 
this work (10^^ to 10^^ W/cm^) seem too low for such events to occur. Intensities 4 to 5 
orders of magnitude greater are thought to be necessary for photo-induced (n,Y) reactions to 
occur via optically-energized electrons. 
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[0122] It is thought that both toroidal and axial magnetic fields are generated 
within laser-induced plasmas and that these fields, near the surface of the ablation target, are 
on the order of 0.6 MG for laser intensities around 10^^ W/cm^. The time development of 
the toroidal field precedes that of the axial field and in so doing provides a trans-radial 
impulse to the ion velocities. The presence of an axial field provides the conditions for a 
plasma centrifuge mechanism to operate in producing the isotope separation we observe. 
The intensity of both fields decreases monotonically away from the surface of the ablating 
target. Data demonstrates that an inverse power law decrease in the axial field occurs over 
distances of about 150 jim from the ablation surface. By extrapolating an initial 0.5 MG 
field to large distances, one can estimate the axial field strength along the path toward a 
deposition substrate on which the isotopically-enriched films are grown. Using a -0.33 
power law dependence for the axial field strength as a function of distance for a laser 
intensity of 2 x lO"* W/cm^, there is obtained a spatially-averaged field strength of 44 kG 
across the 7 cm distance from our ablation target to the deposition substrate. Since the axial 
field is considered to be produced by a dynamo action within the plasma itself, one might 
expect that the field is retained within the plasma plume as it moves toward the deposition 
site. This would provide the necessary geometry for a plasma centrifuge-like condition to be 
established. Such a process provides a mechanism by which enrichment of the light isotope 
occurs on the zero radial spot of the plume. 

[0123] A basic plasma centrifuge equation yields a separated enrichment ratio in 
terms of the radial distance from the center of an axial magnetic field: 
R(r)/Ro=exp(Amo)V/2kT), where R(r) is the ratio of the heavier to lighter ion species at 
radius r normalized to the ratio on axis, Rq. The enrichment depends on the isotope mass 
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difference Am, the angular plasma rotation rate co, Boltzmann's constant k, and the plasma 
temperature T. Based on the data shown in Figure 14, for a factor of 2 enrichment and a 
radial separation of enriched isotopes of 2.5 cm, there is a rotation rate of 3.3 x 10^ rad/sec. 
The value of the equilibrium plasma temperature is taken as 0.5 eV and is based on typical 
electron temperatures measured for similar pulses in titanium using optical emission line 
ratios from spectroscopic emission data. Comparing the present results with the operating 
parameters of conventional plasma centrifuges allows for an estimate of the magnetic fields 
involved in the present experiments. Using this approach, the present minimum plasma 
rotation rate corresponds to a uniform axial magnetic field distributed over the space 
between the ablation target and film substrate of 4 kG. As stated above, the extrapolated 
spatially-averaged field across this region is 44 kG. Thus, the predicted average axial field is 
more than sufficient to cause the isotope separation observed in the present experiments. 
Accordingly, a plasma centrifuge mechanism is consistent with the results presented here. 

[0124] The present invention also provides a method for measuring internal 
magnetic fields of laser ablation plasmas. This is a uniquely valuable scientific result, in 
addition to the use of the inventive technique as a simple and direct means to obtain 
isotopically-enriched materials. It is also useable to obtain other fundamental information 
about laser plasmas derivable from such observations. This includes information about ion 
and electron densities, internal atomic and ionic collisions, ion and electron temperatures, 
ion accelerating fields, and other internal electromagnetic field conditions. From a practical 
standpoint, the direct deposition of engineered isotopically-enriched thin films is an 
important advantage provided by the present invention. 
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Enhanced Ionic Yield 

[0125] Figure 19 contains two graphs showing ionic (nO and neutral (nn) 
concentrations in a laser ablation plume formed by 100 fs pulses in Figure 19A and 10ns 
pulses in Figure 19B, at intensities shown in the data box. As stated in the background 
section of this application, it can be seen in Figure 19 that the ratio of ions to neutrals in a 
single ablation plasma plume is relatively small. Typically that ratio is on the order of 1% or 
less. To achieve a practical method for using this phenomenon as a way to harvest enriched 
isotopes, it would be desired to increase this ionic component to much greater levels. 

[0126] A single laser ablation plume is the result of a nimiber of relatively 
complex processes occurring during the absorption of a laser pulse by a material and 
subsequent expansion of the gaseous/plasma plume. One part of this process involves the 
initial production of a very large percentage of ions if not complete ion content in the vapor 
phase plasma. In other words, the system starts out as a fully ionized electron-ion plasma. 
Because of the high initial density and ensuing collisions in at the early vapor/plasma stage, 
however, a great deal of electronic neutralization occurs as the plasma expands away from 
the surface of the ablating material. This highly neutralized plasma then moves through the 
region of the self-generated magnetic field and the ionic component undergoes the isotope 
enrichment. 

[0127] A plasma having a density which exceeds a certain level, known as the 
"critical density", begins to reflect optical radiation rather than absorb it. The absorption 
properties of a plasma normally increase with the density of free electrons in it. However, 
when the density reaches the critical density the plasma becomes nearly totally reflective. 
This limits the amount of internal energy that can be added to a plasma by a laser pulse. In 
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most cases where an ablation plasma is formed by a laser this critical density is achieved in 
the early part of the laser pulse impinging on the surface of the material. 

[01 28] Referring to Figure 20, another embodiment of the present invention is 
illustrated. In this embodiment, the ionic fraction is increased by first forming the plasma by 
striking the target with an initial laser pulse and then allowing the dense plasma to expand. 
The laser pulse is selected to have an intensity and wavelength sufficient to generate a 
plasma comprising ionized isotopic species and to generate an internal electromagnetic 
accelerating field within the plasma causing spatial separation of the ionized isotopic 
species. 

[01 29] The expanding plasma increases its volume thereby decreasing the plasma 
electron density. This allows the density of the plasma to reduce and to become equal to or 
less than the critical density. During the expansion process a certain degree of neutralization 
occurs in the plasma. This additionally reduces the free electron concentration. 

[0130] After a pre-selected period of time elapses from the time of the initial 
laser pulse, that is, after the plasma is allowed to expand a desired amount, a second laser 
pulse is directed to strike the plasma. The time delay of the second pulse is measured 
relative to the initial laser pulse but is coordinated with the state of the expanding plasma. 
By striking the plasma with a correctly delayed second pulse, the plasma increases its 
absorption of additional energy and its ionic component. 

[01 31 ] The correct time separation between the initial pulse and the second pulse 
is a very important aspect of the present method. If the plasma expansion is allowed to 
progress too long, the plasma density becomes too low to efficiently absorb energy. If the 
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plasma expansion is not allowed to progress long enough, the plasn[ia density does not 
reduce sufficiently relative to the critical density to allow efficient absorption of energy. 

[0132] The present method also includes a technique for identifying the time 
delay desired to maximize the absorption of the second pulse and to thereby improve the 
isotopic enrichment effect. Referring to Figure 21, plasmas of silicon are studied and the 
yield and average energy of the ionic component in the ablation plume as a function of time 
delay between the initial or originating laser pulse and the second or pumping laser pulse is 
shown. The energy and associated flux density in the ablation plasma are a function of the 
percentage of energy extracted from the incoming laser pulse at the target. In Figure 21, the 
central axis of the expanding plume is being examined. 

[0133] Figure 21 shows the results of measuring total ion yield and average 
energy in a pre-formed ablation plume generated by a 120 femtosecond (780 nm) laser pulse 
followed by an identical time delayed pulse. The single pulse at zero delay has an energy 
fluence of 2.2 kJ/cm^ on a beam spot diameter of 42 microns. The two double pulses have 
an energy fluence of 1.1 kJ/cm^ each. The expansion distance is based on a measured 
average ion velocity of 1.9 x 10^ cm/s. Time delays between 0 and 40 picoseconds are 
illustrated. 

[0134] It is seen in Figure 21 that for the conditions of this particular 
experiment, the maximum enhancement above the single laser pulse result occurs at a time 
delay of between 1 and 10 picoseconds, and more particularly between 3.5 and 6.5 

picoseconds, and even more particularly, at approximately 5 picoseconds. This increases the 
ionic yield by a factor of about 2.4 over a single pulse with twice the energy of each of the 
double pulses. In this example it is assvuned that the plasma has expanded from a density 

46 



Attorney Docket No. 2115D-002245 

greater than the critical density (an overdense stage) to the critical density in the intervening 
5 picoseconds. At this density or below, the plasma again becomes an efficient optical 
absorber. 

[01 35] From the foregoing, it can be concluded that the onset of expansion in the 
plasma occurs at less than 0.50 + or - O.lSps. This is because the incident laser pulse 
produces a dense electron-hole solid-state plasma in the time interval during the first laser 
pulse delivery. As a result of the high energy density in this system, a very strong, 
temperature dependent, electron-lattice shock coupling dominates and drives the process of 
converting the solid to a vapor plasma in sub-picosecond time intervals. Subsequent 
expansion of this plasma takes the plasma from an overdense state to a critically dense state 
and then to an underdense state. Critical density occurs when the laser frequency and the 
plasma frequency are equal. 

[01 36] Figure 22 shows an experiment similar to that of Figure 21 except the 
experiment of Figure 22 is for boron isotope enrichment in a pre-formed plasma. In Figure 
22, the isotope enrichment on axis for BN ablation was initiated by a single 120-fs pulse 
(780 nm at 10 Hz) having an energy fluence of 2.2 kJ/cm^ and compared to those initiated 
and pumped by two pulses of half that energy fluence (1.1 kJ/cm^)each and separated by 10 
ps. The total laser intensity was 2x10^^ W/cm^ and the natural abundance is normalized to 1. 
In Figure 22 the normalized ratio of to B^^ is shown for the two cases, as a function of 
charge state. While the stated energy fluences were employed, energy fluences in a range of 
10^ to 10'^ wW are preferred depending on the material selected for the target. 

[0137] It can be seen in Figure 22 that the time delayed, double laser pulse 
method shows an enhancement of enrichment for all charge states. A time delay, in this case 
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10 picoseconds, results in an increase of the total ionic isotope enrichment by about 28% 
over that obtained with a single pulse. For a charge state of +2, the improvement in 
enrichment is about 60%. 

[01 38] The data illustrated in Figures 21 and 22 confirm the concepts as stated 
above relating to the use of time delayed laser pulses for pumping the ionic component of a 
laser plasma to increase isotopic enrichment on the central axis of the plume and further to 
use laser pumping to increase the efficiency of isotope separation on the central axis of an 
ablation plume. 

[0139] Referring now to Figure 23, a second embodiment method for increasing 
the ionic yield even further and achieving even higher isotopic enrichment is described. In 
this embodiment, a plurality of sequentially time delayed laser pulses are directed at the 

plasma after the initial laser pulse. That is, instead of using one time delayed second or 
pumping laser pulse as described in Figure 20, multiple time delayed pumping laser pulses 
are employed. 

[0140] The second embodiment method of Figure 23 is preferred when the 

second laser pulse which pumps the original plasma causes the plasma to again achieve or 
exceed critical density. A third pulse is then directed at the plasma, appropriately time 
delayed, to once again pump the plasma. 

[0141] If the ionic yield is increased by a factor of 2 with each pump pulse and 
the single pulse produces a 1% ionic yield, then a 100% ionic plasma can be reahzed by 
pumping it n times where: n = Lnl00/Ln2 = 6.64. Thus, by pumping the plasma between 6 
and 7 times a 100% ionized condition could be achieved. 
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[0142] If each time delay is on the order of 5 picoseconds, a 100% ionized 

condition can be accomplished in less than 40 picoseconds. With the plasma velocities 
being typically 3 x 10^ cm/sec, the total expansion distance over the 40 picosecond process 
time amounts to 1.4 microns. This distance is much less than the hundreds of microns over 
which the plasma magnetic field extends. As such, abundant time and distance exists for the 
magnetic separation process and associated isotopic enrichment to occur. Even if the plasma 
velocities became as high as 3 x 10^ cm/sec, the total expansion distance will not exceed 14 
microns, thus continuing to be well within the distance needed to interact with the magnetic 
fields. 

[0143] It should also be appreciated that in a multiple pulse operation, the time 
delay between sequential pulses may be varied to yield enhanced results. In other words, it 

may be beneficial to delay each pulse differently. It is also possible that varying each pulse 
duration during the multiple pulse operation could yield enhanced results. As such, the 
skilled practitioner should recognized that by varying individual pulse separation delays and 
individual pulse durations, enhanced control over the operation may be realized. 

Chemical Separation 

[01 44] The above discussion concerns spatially separating isotopes by mass. An 
extension of this principal is chemical separation. That is, different chemicals may be 
separated by mass as well. However, differences in chemistry between the unlike particles 
may alter the separation character efficiency. Isotopes are less complex to separate in that 
they have identical chemical properties. 
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[0145] Referring to Figure 24, a graph is provided illustrating elemental 

enrichment of nickel (Ni) and copper (Cu) as a function of angle in the ablation plume. The 
total enrichment ratio is shown first for each angle, followed by each of the ten charge states. 
The circles and solid line indicate the total enrichment ratio. The actual yield of each charge 
state and angle is not shown in the graph. It should be noted that the sample ratio of 45% Ni 
to 55% Cu was normalized to 1. 

[0146] Nickel and copper are both metals which are next to each other on the 
periodic table. Nickel and copper in that regard are both metals with similar chemical 
properties, and thus separate in a similar way to that of isotopes. Figure 24 shows that the 
lighter element (nickel) is more abundant at 0 degrees, just like the case in isotope 
separation. 

[0147] Figure 25 is a graph illustrating elemental enrichment of strontium (Sr) 
and titanium (Ti). Strontium's atomic mass is 88. Titanium's atomic mass is 48. The 
chemistry of strontium and titanium is quite different in the sense that they have very 
different electronic structure leading to very different ionic yields and velocity distributions. 
Thus, they do not separate in the same way as isotopes. In Figure 25, it can be seen that the 
heavier element (strontium) is actually more abundant at 0 degrees, in contrast to the isotope 
case (where the lighter isotope is more abundant at 0 degrees). 

[01 48] The method for conducting chemical separation is the same as for isotope 
separation. That is, the chemicals are separated based on mass separation in a laser ablation 
plume presumably by self-generated electric and magnetic fields. The only difference is the 
target. For Nickel/Copper, a Ni;Cu alloy target was used (45% Ni, 55% Cu). For 
Strontium/Titanium, a SrTi03 target was used (strontium titanate). Even though the method 
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is the same, because the two particles in Figure 25 have very different chemical structures, 
the separation character efficiency is different. In the case illustrated in Figure 25 (Sr/Ti), 
the efficiency is inverted compared to that of isotopes. Other chemical pairs will have varied 
efficiencies and may differ from that of Figure 24 and 25. 

[0149] As stated above, the method and system apparatus used for chemical 
separation is essentially the same as the method and system apparatus used for isotope 
separation. As such, the method includes directing a laser beam onto a surface of a target. 
The laser beam has an intensity and wavelength sufficient to generate a plasma comprising 
chemical elements and to generate an internal electromagnetic scattering field within the 
plasma. This causes spatial separation of the chemical elements. 

[0150] If desired, the spatially separated chemical may be deposited on a 
substrate. The deposit would be formed with a first region which is relatively rich in a 
selected chemical element. A second region would be relatively lean in the selected 
chemical element. 

[0151] In the above method, the intensity may be in a range of about lO' 

watts/cm^ to lO'^ watts/cm^. The wavelength may be in a range of about 200 nanometers to 
about 10 micrometers. Further, the laser beam may comprise one or more pulses each 
having a duration in a range of nanoseconds to femtoseconds. 

[01 52] Also if desired, an additional laser beam may be directed onto a selected 
one of the first and second regions at an intensity and wavelength sufficient to generate a 
second plasma comprising chemical elements and to generate an internal magnetic field 
within the second plasma. This causes spatial separation of the chemical elements of the 
selected region. Thereafter, the spatially separated chemical elements may be deposited on a 
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substrate, whereby a second deposit is formed having a third region with a chemical 
distribution different from the selected region. This process may be repeated indefinitely to 
produce a highly enriched separation. 

Enhanced Ionic Characteristics 

[0153] Femtosecond laser-pulse absorption above and below the critical density 

regime in ablation plasmas by means of a pair of identical collinear laser pulses separated by 
a few picoseconds is now discussed. The second laser-pulse modifies the ionic 
characteristics of the pre-formed plasma, such as ion yield, ion energy and average charge 
state. Maximum effects occur when the density of the initial plasma drops to just below 
critical density. The enhancement of ion yield is significant compared with that produced by 
a single pulse with twice the energy of each individual double pulse. 

[0154] Application of ultrafast lasers to materials synthesis and processing offers 
interesting new opportunities as compact ultrafast laser systems become more affordable and 
reliable. Particularly, ultrafast, micromachining, optical emission spectroscopy and pulsed- 
laser-deposition (PLD) are examples of potential benefits and interests since one can take 
advantage of the adjustable properties of ablation plumes. These plumes, and the associated 
plasmas, can be pre-tuned by varying laser pulse parameters such as pulse width, intensity 
and contrast. 

[01 55] Interacting the initial plasma with a time-delayed second laser pulse can 

further modify the characteristics of an ablation plume. For example, an experimental 
investigation was conducted into the enhancement of ionic properties in the plasma portion 
of ablation plumes from a silicon target resulting from two identical 110 fs laser pulses 
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separated by a few picoseconds. Time-resolved measurements were made of ion flux, ion 
energy and average charge state to verify that the ionization level and electron density 
condition of a plasma, initially formed by a single pulse, can be subsequently modified by 
the pumping action of a secondary time-delayed pulse and that critical density effects could 
be usefully manipulated. 

[0156] The laser source used in this experiment was a Ti:Sapphire CPA laser 
system operated at lOHz. The laser provided 110 fs, 36 mJ laser pulses at a wavelength of 
780nm. For this experiment a single laser pulse was split into two identical pulses with a 
variable time delay using a Michelson interferometer. Resolution of time delay was 0.2 ps. 
Both pulses were p-polarized and focused collinearly on a silicon target using a 750mm focal 
length lens. The incident angle was 45 degrees. The ablation target was mounted inside a 
vacuum chamber with starting pressure of 5x10"^ Torr. 

[0157] In this experimental setup, time-of -flight (TOP) ion signals can be 
measured using a curved sector electrostatic energy analyzer (with multi-channel plate, 
MCP, detector) to obtain energy distributions as a function of charge state. Alternatively, the 
ion flux can be allowed to proceed undeflected into a similar MCP on a straight through port 
of the analyzer. In the later case an undifferentiated ion flux density is obtained. The 
entrance of this detector system is normally located at a distance of 150 cm from the ablation 
target in a separate, differentially pumped, vacuum chamber. 

[0158] Total ion flux densities, from the straight through port, were obtained by 
time integration of the TOP signals. Energy distributions, and average ion energy, were 
obtained from the straight through port by converting the TOP data to kinetic energy. Ion 
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yields as a function of charge- state were obtained using the curved- sector portion of the 

energy analyzer. 

[0159] Figure 26 shows the total area density of ions in the ablation plume 
produced by a pair of pulses, as a function of time delay between the pulses, and the laser 
pulse intensity. The data points at zero delay are the result of a single pulse with twice the 
energy as that contained in each of the individual double pulses. It also represents the result 
of delivering the double pulses with zero time delay. 

[01 60] Two effects are observed in the data of Figure 26a when the second pulse 
maximally interacts with the plasma formed by the initial pulse. First, the ion yield is three 
times higher than that produced by a single pulse. This is a result of further ionization of the 
plasma by the second pulse. Second, ion average energy is increased as a result of the second 
pulse. Ion kinetic energy is 9keV using the double pulses separated by 5 ps compared with 
6keV using a single pulse. It is probable that these results are a consequence of the bipolar 
electrostatic field, responsible for ion acceleration, being sustained and reinforced by 
additional electrons from the second laser pulse so that ions are experiencing a stronger field 
for a longer period of time. 

[0161] In Figure 26(a) the maximum enhancement occurs at a 5 ps delay 
between the initial and subsequent pulse. This impUes that 5 picoseconds are desired to 
allow time for sufficient expansion of the initially over-dense plasma, produced by 
IxlO^^W/cm^, to develop and drop below critical density. Once the plasma density reaches 
critical density, and goes slightly below, it can again become an efficient optical absorber. 

[01 62] Critical density is defined here in the usual way as, ric = mco^/47ie^ where 
nc is the electron density, m the electron mass, co the laser wavelength, and e the electronic 
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charge. For the wavelength of our experiment (780 nm), the critical density is 1.7 x 10^^ 
electrons/cm^. The plasma functions as a highly reflective media above this density and is 
absorbing just below. At a sufficiently low density the plasma becomes essentially 
transmissive, 

[0163] The absorption properties below critical density increase with electron 
density up to the critical point whereupon the plasma becomes reflective and absorbs 
inefficiently in a shallow skin depth. Based on these arguments and a self-similar expansion 
model for the ablation plasma an estimate can be made of the initial plasma density prior to 
expansion. For the plasma of Figure 26a to achieve critical density in 5 ps, with a known 
average velocity of 1.8 x 10^ cm/sec, the plasmas initial density would have to be 6 xlO^^cm" 
^. This is approximately 3.5 times higher than critical density, at the wavelength being used. 

[0164] The onset of expansion is seen to occur in Figure 26 at less than 0.5 ps. 
This implies a very strong coupling of laser energy to the plasma in sub-picosecond time 
intervals. It can be inferred from this, that as a result of the extremely high energy-density in 
the electronic sub-component of the plasma, a strong electron temperature-dependent shock 
coupling to the lattice occurs. This process results in the conversion of the solid material to 
a vapor-state plasma in sub-picosecond time intervals. 

[01 65] Ion yields as a function of pulse separation are plotted in Figures 26 (b) 
(c) and (d) using 75%, 48% and 19% of maximum laser intensity respectively. In comparing 
these data, the ion yield and average energy enhancement factors systematically decrease 
with laser intensity. This is presumably because the initial plasma densities for the higher 
intensities are significantly above critical density and, after sufficient expansion time, can be 
pumped back up to the higher ionic density and energy level. 
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[01 66] However, those at the lower intensities are closer to the critical level and 

present a more tenuous plasma to the secondary pulse. In these cases, especially for Figure 
26(d), which represents results produced by laser pulses with 19% of the maximum 
intensity, the initial plasma densities are probably close or below critical density. Under 
such circumstances, as the plasma expands the density decreases further below critical 
density and therefore the second laser pulse, in penetrating such a dilute plasma, has limited 
effects on both ion yield and ion kinetic energy. 

[0167] In further support of this argument, it is clearly seen, by the diagonal 
dashed line running from Figure 26(a) through Figure 26(d), that less time delay is required 
to reach the maximum enhancement as laser intensity decreases. This peak shifting effect is 
considered to be a result of the initial density of the laser-produced plasmas being less at the 
lower laser intensities and ultimately being sufficiently low so as to approach the critical 
density. 

[0168] If an approximately linear decrease of initial plasma density with laser 
intensity is assumed, the initial 6 x 10^^ cm^ density of Figure 26a is reduced to 1.2 x 10^' 
cm^ at the 19% level in Figure 26d. As pointed out above, this is slightly below the critical 
density of 1.7 x 10^^ cm'^. The desired density therefore is between about 65% critical 
density to about 100% critical density. In the case of Figure 26d, the second pulse 
essentially penetrates the initially formed plasma at all delays and adds to the single plume 
density to provide a relatively constant ion flux density at all time delays. 

[0169] Figure 27 shows the distribution of charge states for the single and 
double pulses at maximum laser intensity. The average ion charge in the plume for the case 
of a single pulse is 2.5 and it increases to 3.1 in the case of a double pulse with a 5ps time 
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delay, and 3.0 for a lOps delay. The inset in Figure 27 shows the general single mode shape 
of the energy distribution for the higher intensities used in Figures 26 and 27. 

[01 701 It can be seen in Figure 27 that absorption of the second pulse results in a 
significant increase of the higher charge state ions (+5 to +7) while having modest or no 
effect on the number of low charge state ions. The collisional-radiative equilibrium (CRE) 
plasma model can be used to describe such plasma ionization dynamics. This model 
normally shows an increase in the number of higher charge-state ions at the expense of lower 
charge states as the plasma temperature increases. This is not the case for the data of Figure 
27. 

[01 71 ] The increase of average ionic charge does not occur at the expense of the 
low charge state ions but rather by only generating a larger number of high charge state ions. 
These additional contributions to the high charge states may be associated with the 
absorption of the second pulse predominantly in the leading edge of the expanding plasma, 
which is where the high-velocity high-charge state ions normally reside. One can expect 
resonance absorption to be in effect with p-polarized light, which is incident at a 45-degree 
angle, resulting in preferential absorption near the front critical surface of the expanding 
plasma. 

[01 72] If the laser intensity is further reduced below that of Figure 27 a new and 
more subtle effect emerges. Figure 28 shows ion energy and ion yield in an ablation plume 
produced by laser pulses with 6% of maximum intensity (6x10^"^ W/cm^) as a function of 
time delay. A characteristic ion energy distribution is shown in the inset of Figure 28 and 
consists of two components rather than the single large component shown in Figure 27. As 
shown before, the low energy component, centered at lOOeV, is relevant to the 
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hydrodynamic expansion, while the high-energy component, centered at IkeV, is driven by 
the internal bi-polar plasma electrostatic field. 

[01 73] As shown in both Figures 27(a) and (b), the overall (total) ion energy and 
yield are essentially unchanged when using single pulse or double pulse. However, when the 
two components (high energy and low energy) are considered separately, a significant time 
delayed signal enhancement is found for the lower-yield low energy component. 

[0174] These results suggest that the density of the hydrodynamic component in 
the plasma produced with the 6% laser intensity is closer to (or above) critical density 
compared to that of the high-energy component. The time-delayed second laser pulses are 
therefore able to more effectively modify characteristics of the low energy ions. 
Furthermore, since the low energy component expands relatively slow, the range of time 
delay giving ion yield enhancement is wider than that shown in Figure 26(a). 

[01 75] Since high-energy ions are dominant in the yield for the overall (total) ion 
distribution in terms of percentage, the total ion response does not reflect the enhancement 
associated only with the low energy component. The broad energy distribution of the high- 
energy component in Figure 28 suggests that it is more tenuous than the low energy 
component, which has a narrower energy distribution. In these expanding plasmas the high 
velocity ions spread out ahead of the low energy ones to produce a spatially distributed 
concentration distribution. The single mode spatial distribution of Figure 27 results from 
the merging of the two component bi-modal distribution for plumes generated at those 
higher laser intensities. 

[0176] In summary, double pulses separated by picoseconds effectively increase 
ion yield, ion kinetic energy and average charge in the plasma compared with those produced 
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by a single pulse with twice the energy of each individual pulse. The maximum 

enhancement occurs whenever the plasma reaches a level close to or slightly below critical 
density so that it can effectively absorb the second laser pulse. The time-resolved 
enhancement is therefore an interesting method for performing plasma expansion processes. 
Furthermore, since the internal electromagnetic field can be sustained or reinforced by the 
second pulse, effects such as isotope separation can be enhanced using the double pulse 
technique. 

[0177] An operational approach can be envisioned where multiply sequenced 
pulses are used to maximize the ionic component of the plume, where the number of pulses 
may be on the order of 5 to 10, each one taking the plasma back to critical density in a 
nominal 5 ps time interval. A highly ionized plume could then be achieved in a time interval 
of 30 to 50 ps. 

Angular Distribution of Mass Enrichment for Heavy Ions 

[0178] The following concerns spatial mass enrichment of chemical and isotopic 
ions in laser plasmas generated in ablation plumes and particularly from an alloy of Ni/Cu 
using ultrafast laser pulses. It is shown that ions of chemicals with neighboring masses can 
be enriched in these plasma plumes similar to the way in which isotopes are enriched. 
Comparison of angular distributions for the ratios of the elements nickel (mNi=58.71) and 
copper (mcu=63.54) are made with those found for similar alloy plasmas subjected to 
conventional plasma centrifuge processes. The results for both chemical species and nickel 
isotopes support a plasma centrifuge mechanism for the spatial enrichment of these heavy 
ions in the ablation plasma plumes. 
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[0179] As described above, spatial isotope enrichment for boron can be achieved 
on the central axis zone of an ultrafast laser ablation plume. Similar results have been 
demonstrated by the present inventors for a wide range of heavy ions using a variety of laser 
pulse durations and wavelengths. Although some preliminary data on angular distributions 
for boron enrichment in ablation plasmas was previously collected, its support for a plasma 
centrifuge mechanism was somewhat inconclusive. The following shows that such a 
centrifuge mechanism is rather strongly supported by direct comparison of chemical and 
isotopic enrichment as a function of angle about the normal to the ablation surface. 

[0180] The chemical separation is based on examining plasmas generated from 
alloys of nickel and copper (Ni 45%, Cu 55%) where the average alloy mass is 61.37 mass 
units and the species mass difference is 4.85 mass units. In this case 6ni/m is 4.85/61.37 = 
7.9 X 10"^. For the case of boron isotope enrichment the 5m/m was 1/10.81 = 9.9 x 10"^. 
Since these numbers are similar, it might be expected that the alloy system would show a 
similar behavior to that of the boron isotopes, where now however two different chemical 
species are studied rather than two different isotopes of the same chemical. 

[0181] Going beyond the similarity in differential mass ratio, it is instructive to 
examine the isotopes of nickel (^^Ni and ^°Ni), which are representative of a relatively heavy 
element. In this case 5m/m is 2/58.98 = 3.4 x 10'^ which is about a factor of 3 smaller than 
for boron isotopes. It is instructive to examine this case since if a plasma centrifuge process 
is in effect, the relative spatial mass separation for nickel isotopes should be similar to that 
of boron, since a centrifuge mechanism is inherently independent of 5m/m ratios. Most 
other mass separation mechanisms, and especially hydrodynamic mechanisms, tend to 
rapidly become ineffective as this ratio decreases. 
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[0182] The angular distribution of enrichment about the normal to the ablating 

surface is an important physical property to be examined. If a centrifuge process is in effect, 
this distribution should follow the centrifuge equation in relation to angle and be 
conformally related to the results obtained from a conventional cylindrical centrifuge in 
relation to radius. The transformation between these two coordinate systems is shown in 
Figure 29 where the model of the present invention is depicted. 

[01 83] As described above, the generation of a laser-plasma by an ultrafast laser 
pulse involves the absorption of laser energy by the electrons in the solid and the transfer of 
that energy to the atoms and ions in the material structure. This happens on a very short time 
scale on the order of a picosecond or so. The expansion of the ablation plasma occurs 
initially as a one-dimensional expansion process moving along the z-axis normal to the 
surface of the ablating material. The plasma converts to a three dimensional expansion 
process after a distance that is on the order of the laser beam spot dimension. 

[0184] The isotopic and/or chemical mass separation effects are expected to take 
place in this one-dimensional zone. This is where the magnetic fields are the strongest and 
where centrifugal action can occur as a result of high-velocity ionic interaction with the 
magnetic fields. In fact, it can be shown that the plasma ions are able, in principle, to rotate 
at the cyclotron frequency in the magnetic fields of the ablation plasma where these field 
strengths are on the order of 100 K-gauss to 1 M-gauss. 

[01 85] The equilibrium distribution of the ions in the one-dimensional cylindrical 
zone are similar to that achieved in a conventional plasma centrifuge. When the ions move 
into the three-dimensional expansion zone they move out of cylindrical symmetry and into 
angular symmetry which is characteristic of a pulsed laser plasma. 
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[0186] Figure 29 illustrates a coordinate transformation showing the relationship 

between the one-dimensional cylindrical zone of the plasma centrifuge and the three- 
dimensional expansion zone with the angular distribution associated with the ablation 
plasma. Figure 29 shows an idealized behavior for how these ions would undergo this 
geometric transformation process. Those ions on the central axis of the one-dimensional 
cylindrical zone remain on the central axis of the three-dimensional conical expansion zone, 
whereas the ions nearer the radial perimeter of the one-dimensional cylindrical zone expand 
out to the largest angles of the three-dimensional conical expansion zone. 

[0187] The generalized centrifuge equation, valid for isotopic separation in 
mechanical gas centrifuges or plasma rotating devices, is given as R(r)/Ro = exp(- 
AmcoV/2kT) where R(r) is the ratio of the lighter to heavier ion species at radius r 
normalized to the ratio on axis, Rq. The other terms in the equation are Am, the mass 
difference of the species being separated; co, the angular rotation rate of the centrifuge 
process; r the radius of the cylindrical zone where enrichment is being observed; k, 
Boltzmann's constant, and T the absolute temperature of the rotating gas or plasma. In the 
case of a magnetically driven plasma centrifuge, the rotation rate can be one of two distinct 
angular frequencies, one at the cyclotron frequency and one at a specific lower rate, 
proportional to the E x B drift rotation. Derivation and discussion of these two rigid-rotor 
rates is known in the art. 

[01 88] It is instructive to note that the enrichment process is independent of the 
actual mass involved and depends only on the difference in mass of the species to be 
separated. In the experimental arrangement used herein, the mass separation is expected to 
occur in the one-dimension cylindrical expansion zone. These effects are carried forward 

62 



Attorney Docket No. 2115D-002245 

into the three-dimensional conical expansion zone. Figure 29 shows a schematic 
representation of this process. 

[0189] Mass separation, by the centrifuge process, results in an 
ordering/distribution of the high and low mass species along the radius of the one- 
dimensional cylindrical expansion zone. This ordering of rotational energy is balanced by the 
randomizing thermal energy kT to result in a steady state distribution. 

[01 90] In an ordinary mechanical or plasma centrifuge device, the radial ordering 
is observed directly at the end of the cylinder. In the illustrated embodiment, the plasma 
transitions from the one-dimensional cylindrical expansion zone to the three dimensional 
conical expansion zone a short distance from the ablation target. It is therefore necessary to 
map the one-dimensional radial ordering in r to a three dimensional angular distribution in 
0. 

[0191] For this purpose, an elementary conformal mapping function is used to 
convert the coordinate r to the angular dimension rni(0 / 0 m) where rm is the maximum 
radius in the one dimensional zone, and is equal to the radius of the laser beam spot, and 0 m 
is the maximum angle to which rm maps. This maximum angle is determined from 
experimental observations. Figure 29 demonstrates how this transformation is effected. The 
angular expression for the centrifuge equation is discussed further below. Experimental 
angular distribution data for the Ni/Cu system, based on these concepts, are presented below. 
In addition, a discussion of the parameters in the angular form of the centrifuge equation is 
also presented as support for the magnetic centrifuge model operating in the one-dimensional 
expansion zone and manifesting itself in the three dimensional expansion data. 
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[0192] Figure 30 shows an ion energy density distribution for ^^Ni ions as a 
function of charge state. These results are similar to those found with boron as described 
above. The pattern in Figure 30 is very characteristic of energy distributions in plasmas 
formed by ultrafast laser ablation plumes. The lowest charge state ions have the highest 
yield and lowest energies, whereas the highest charge state ions have the highest energy and 
lowest yields. These data are obtained from an ion energy analyzer operated in a time of 
flight mode. By collecting data in this fashion it is possible to extract isotopic or mass ratios 
as a function of ionic charge state. Examining the data in that fashion helps to understand 
details of the model and provides additional insight to the enrichment process 

[0193] As described above, Figure 25 illustrates the angular distribution of the 
ionic ratios representing the chemical mass enrichment in the ablation plasma from the alloy 
of Ni/Cu. The alloy ratio of 45% Ni and 55% Cu is normalized to unity and the angular scan 
represents either enrichment or depletion of the Ni ions relative to that of Cu. It is seen that 
at zero angle the ratio is increased by about a factor of 2.5 above the alloy composition and 
at 60 degrees the plasma is depleted by an approximately equivalent amount. The behavior 
illustrated in Figure 25 is consistent with what would be expected from a centrifuge process 
as depicted in Figure 29. 

[0194] A direct comparison can be made between Figure 25 and similar 
experiments done with a conventional plasma centrifuge. In so doing, it is observed that the 
enrichment using ultrafast ablation is about 2 Va times greater than that observed with a 
conventional magnetic centrifuge. This presumably is a consequence of the much higher 
magnetic fields available in these laser plasmas. 
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[0195] Figure 31 illustrates the angular distribution for the isotope ratio of 58Ni 

to 60Ni as a function of ionic charge state. The general pattern of light isotope enrichment 
on axis is seen in this data as well as the decrease of this enrichment with increasing angle. 
The solid circles and solid black line in the Figure 31 shows the weighted average over all 
charge states for the enrichment.at each of the angles measured, which were 0, 15, 30, and 60 
degrees. 

[0196] The isotopic ratio is seen to come back to slightly less than natural 
abundance at 60 degrees. The depleted ^^Ni ratios are presumably at the higher angles and 
compensate for the higher values at zero degrees in a fashion similar to that of Figure 25. 

[0197] An unusually high enrichment of the charge state 8 and 9 ions are seen in 
Figure 31. These ions represent a small amount of the total ionic yield and therefore do not 
contribute heavily to the total ionic enrichment. However, the presence of these ions 
suggests that some interesting physics is occurring to provide enhanced enrichment for these 
particular species. A similar effect was seen with boron, where charge states 3 and 4 had the 
highest light isotope enrichment but a low yield. 

[0198] Figure 32 illustrates the angular distribution of ionic yield as a function 
of charge state for the combined 58Ni and 60Ni ions. It is seen that the highest yield results 
from charge states 1 through 3 with a maximum near the central axis and a minimum at 60 
degrees. The remainder of the higher charge state species exhibit much lower total yield at 
all angles. 

[0199] Figures 25 and 31 provide a set of data that can be used to evaluate the 
magnetic centrifuge model in relation to the model diagram of Figure 29. It is shown here 
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that a highly self consistent set of numerical values concerning rotation rates, magnetic field 
strength, and ion velocities emerge from the model and the data presented. 

[0200] Starting with the centrifuge equation, as transformed for angular 
expansion, R(0)/Ro = exp(-5mco^rin^0^/0ni\T) where all terms are defined as stated above. 
Inverting the centrifuge equation and solving for the rotation frequency gives: co = [- 
lM{Rie)fRo}@j2kTf6mrJ&Y^- Using values from Figure 25; Ro = 2.2, R(30°) = 0.80, 
0in = 60^, and 0 = 30°. Taking 5m = 4.85 mass units for the Ni/Cu system and r^ = 150 
microns (beam spot radius), the angular rotation rate becomes &> = 7.86 x lO^T^'^ where T is 
the plasma temperature in absolute degrees. 

[0201] For the centrifuge mechanism to produce an ordering in the rotational 
energies of the plasma ions and therefore to achieve an appropriate spatial separation, the 
angular rotation rate must be large enough to balance against the randomizing motion 
induced by the mean temperature of the plasma. The angular rotation rate scales as the 
square root of the plasma temperature and must, if the separation is to continue, be driven 
higher by an appropriate magnetic field as the temperature of the plasma increases. This is 
an important concept in determining details of the plasma mechanism involved in producing 
the desired results. 

[0202] The initial temperature of a plasma, formed by pulses such as those 
described herein for metallic samples is on the order of 200 ev. After about 20 ns of 
expansion time this temperamre drops to about 1 ev. In the early stages of plasma formation 
then, the temperature is on the order of 2.3 x 10^ K requiring a plasma rotation rate of 1.2 x 
10^ rad/sec to achieve the results of Figure 25. After 20 ns of expansion this rate can be as 
low as 7.5 x 10^ rad/sec. 
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[0203] It is recognized here that the centrifuge equation represents a steady state 
plasma that executes its rotation in a rigid-rotor motion. However, for the sake of 
discussion, and to establish some order of magnitude sense of the physics involved, single 
particle motion is examined for clues as to how the process works. 

[0204] Using the basic equation for the cyclotron rotation rate coc of a charged 
particle q with mass m, moving orthogonal to a magnetic field B, ooc = 27rqB/w, one finds, 
for an ion of charge state 1, a rotation rate of 9.7 x 10^ rad/ sec per megagauss (Mgauss) of 
magnetic field strength. Comparing this to the required rotation rate at the beginning of 
plasma formation implies the presence of a longitudinal magnetic field on the order of 1.2 
Mgauss to initially drive the centrifuge mechanism. Field strengths of this magnitude are 
known to exist in ultrafast laser plasmas. 

[0205] As the plasma expands forward, the temperature drops and the 
corresponding field strength drops. Public data suggest that the longitudinal magnetic field 
drops in an exponential like fashion with a 1/e decrease to 0.44 Mgauss occurring over an 
approximately 100 micron expansion distance. After 20 ns the required rotation rate implies 
the presence of a magnetic field on the order of 70 Kgauss. The plasma will have begun to 
move into the three-dimensional expansion zone by this time. 

[0206] To further explore the concepts of single particle motion, Larmor orbital 
radii for cyclotron motion are given as rc =v/(Dc where v is the orthogonal linear velocity of 
the charged ion. If the rotational motion is confined to occur within the laser beam diameter, 
corresponding to the one-dimensional expansion zone, then the linear orthogonal velocity 
associated with the experimental rotation rates of 1.2 x 10^ rad/sec is 2.9 x 10^ cm/sec, 
corresponding to an average ion kinetic energy of 0.25 Kev for th elements shown in Figure 
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29. It is seen in Figure 29 that this is approximately the mean kinetic energy measured for 
the charge-state one ions, which are the most abundant species, observed in the plasma. It 
should be mentioned however, that the data of Figure 29 are for the longitudinal velocities 
and are therefore only representative of the actual orthogonal velocities in the plasma. There 
is nonetheless good order of magnitude agreement. 

[0207] It needs to be emphasized that care should be exercised concerning too 
great a use of the single particle model to describe the collective plasma motions discussed 
herein. However, a further point of interest is seen in Figure 31 where the isotopes of nickel 
appear to reveal a higher enrichment for the higher charge states, implying higher single 
particle rotation rates for the high charge state ions, as prescribed by the cyclotron equation. 

[0208] In summary, these experiments provide data that support a magnetic 
centrifuge model for the ionic enrichment of light isotopes on the central axis of an ultrafast 
laser ablation plume. The results are similar in relation to the chemical enrichment of nickel 
over copper from a Ni/Cu alloy as well as for the ratios of the isotopes of 58Ni/60Ni from 
the nickel component of the plasma plume. The work further confirms earlier results on 
boron isotopes and on axis data for other heavy ions. The support for a magnetic centrifuge 
model is found in directly comparing ablation plasma results with the spatial enrichment data 
published for conventional wire wound magnetic centrifuge devices. 

Clusters and Aggregates 

[0209] Often, small clusters or aggregates form in a laser ablation plume as the 
plasma expands and cools. These aggregates are not always beneficial and may be a 
detriment to the process. For example, such an aggregate can interfere with the isotope 
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enrichment process by extracting atoms or ions from the plasma and not letting them 
participate in the enrichment process. 

[0210] To prevent such aggregates from forming or to vaporize them once they 
form, a second laser pulse similar in nature to those described above, that is appropriately 
timed, can be used. The time delay required can be determined based on testing of particular 
target material and laser pulse interactions. The second laser pulse may be provided by a 
second laser but is preferably split off from the main laser pulse and delayed in time with, for 
example, a Michaelson interferometer. In addition, multiple sequential pulses may also be 
employed. Tn the case of aggregates, by maximizing the ionized fraction, the aggregates 
m should automatically be minimized or eliminated with the secondary pulses. 
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[0211] While preferred embodiments, forms and arrangements of parts of the 
invention have been described in detail, it will be apparent to those skilled in the art that the 
disclosed embodiments may be modified. Therefore, the foregoing description is to be 
considered exemplary rather than limiting. 
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